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Q ! ABSTRACT 

^ : 

[/3 ' Context. The infrared [Nell] and [Nelll] fine structure lines at 12.81/im and 15.55^m have recently been theoretically predicted to trace the 
. circumstellar disk gas subject to X-ray heating and ionization. 

■ Aims. We observationally investigate the origin of the neon fine structure line emission by comparing observations with models of X-ray 
I irradiated disks and by searching for empirical correlations between the line luminosities and stellar and circumstellar parameters. 
^ Methods. We select a sample of 28 young stellar objects in the p Ophiuchi star formation region for which good quality infrared spectra 
("""^ and X-ray data have been obtained, the former with the Spitzer IRS and the latter with the Deep Rho Ophiuchi XMM-Newton Observation. 

, We measure neon line fluxes and X-ray luminosities; we complement these data with stellar/circumstellar parameters obtained by fitting the 
f — spectral energy distributions of our objects (from optical to millimeter wavelengths) with star/disk/envelope models. 

Results. We detect the [Ne II] and the [Ne III] lines in 10 and I cases, respectively. Line luminosities show no correlation with X-ray emission. 
The luminosity of the [Nell] line for one star, and that of both the [Nell] and [Nelll] lines for a second star, match the predictions of 
f"**) published models of X-ray irradiated disks; for the remaining 8 objects the [Nell] emission is 1-3 dex higher than predicted on the basis of 
their Lx- However, the stellar/circumstellar characteristics assumed in published models do not match those of most of the stars in our sample. 
, Class I objects show significantly stronger [Ne II] lines than Class II and Class III ones. A correlation is moreover found between the [Ne II] 
' line emission and the disk mass accretion rates estimated from the spectral energy distributions. This might point toward a role of accretion- 
generated UV emission in the generation of the line or to other mechanisms related to mass inflows from circumstellar disks and envelopes 
and/or to the associated mass outflows (winds and jets). 
^ Conclusions. The X-ray luminosity is clearly not the only parameter that determines the [Ne II] emission. For more exacting tests of X-ray 
irradiated disk models, these must be computed for the stellar and circumstellar characteristics of the observed objects. Explaining the strong 
[Nell] emission of Class I objects likely requires the inclusion in the models of additional physical components such as the envelope, inflows, 
and outflows. 
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1. Introduction 

The first million years in the formation of a low-mass star are 
characterized by several complex and still not fully understood 
phenomena involving the circumstellar envelope, the disk, and 
the central protostar, e.g., envelope and disk mass accretion, 
outflows, disk evolution including planet formation, star/disk 
magnetic interactions, and other manifestations of the mag- 
netic field such as the intense X-ray emission from hot mag- 
netically confined plasma. The first studies of Young Stellar 
Objects (YSOs) have often neglected important interactions be- 



tween these phenomena. A noteworthy example is the effect 
of X-rays from the central object on the surrounding molec- 
ular cloud, the accretion envelope, and the disk. YSOs, in- 
deed, have stronger X-ray emi ssion than main sequence stars 
jFeigelson & Montmerl3l 19991) . 
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The origin of magnetic phenomena in YSOs, and 
of their X-ray emission in particular, is an intriguing 
and still poorly understood problem. Renewed interest 
in YSO X-ray emission comes from the recent recogni- 
tion that X-rays may ionize and modify in several other 
ways th e thermal and chemical structure of star forming 
clouds dLorenzani & Fallal l200lh . circumstellar disks (e.g. 
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Glassgold et alJll997t lllgner & Nelsoiil l2006t iMeiierink et al 



2008t iGorti & Hollenb ach 200 81: lErcolano et al. 2008), and 
planetary atmospheres ( iCecchi-Pestellini et al.ll2006 ). We here 
focus on the response of YSO disks to X-ray irradiation by the 
central (proto)star. 

The structure and temporal evolution of circumstellar disks 
is of paramount importance for the understanding of the star- 
and planet-formation process. The structure of the gas-phase 
component, by far more massive than the dust component, is 
particularly uncertain being the most affected by high energy 
radiati on (far/extreme ultraviolet and X-ray), loiassgold et al.l 
(Il997h have shown that, in a protostellar disk illuminated by a 
central X-ray source. X-ray ionization dominates over that of 
galactic cosmic rays, giving rise to a vertically layered ioniza- 
tion structure with an outer active surface and a m ostly neu- 
tral inner dead z one (IGlassgold et al.l Il997l: iGamm ie 1996) 



Ilgner & Nelsoiil (1200^ calculated the disk ionization struc 



ture as a function of the X-ray luminosity and emitting plasma 
temperature, and found that the disk is divided into three dis- 
tinct radial zones: an inner active region, a central region where 
the depth of the dead-zone depends on the X-ray spectral and 
temporal characteristics, and an outer region with non-variable 
dead-zone. In addition to ionization. X-rays can heat the ex- 
ternal layers of disk atmospheres, as shown by lGlassgold et al ' 



( 120041) who predict temperatures of the order of 5000 K. 

Theoretical calculations depend critically on several in- 
gredients: the disk model, the chemical network, the spectral 
and temporal characteristics of the X-ray source and its as- 
sumed spatial location with respect to the disk. Observational 
tests are therefore highly desirable and could help constrain 
the model assumptions. The lines of ionized neon are par- 
ticularly useful as a proxy of the effect of high energy ra- 
diation, as its P* and 2"'' ionization potential are 21.56 and 
41 eV, respectively and therefore ionization can occur only by 
photons in the EUV and X-ray range. Moreover, due to its 
clos ed sh ell configuration, the Ne che mistry is particularly sim- 



ple. IGlass gold et al. (2007), Meiieri nk et all (12008, hereafter 



MGN08), IGo rti & Ho llenbach) (120081 hereafter GH08), and 
Ercolano et al.i (.2008.) have recently calculated the strength of 
fine structure emission lines fro m ionized neon o r iginati ng in 
a disk exposed to stellar X-rays. Glassg old et al.l (l2007b esti- 
mate that in low-mass YSOs the ionization of neon is domi- 
nated by X-rays, because the photospheres of these stars emit 
few UV photons and cosmic rays are removed by the strong 
winds. Nell and Nelll ions, predominantly resulting from 
K-shell photoionization of neutral neon by X-rays with en- 
ergy E>0.87KeV, give rise to magnetic dipole transitions at 
12.81/^m and 15.55/im, respectively. The predicted line lumi- 
nosities are, for the reference disk/star models of MGN 08 and 
GH08, and for X-ray luminosities of ~ 2 • lO'^^erg s^^, of 
the order of 10^* and lO^^erg s^^ for [Nell] and [Nelll], re- 
spectively. MGN 08 predict that the line luminosities increase 
with X-ray lum i nosity following a steeper-than-linear relation. 
Ercolano et aP (l2008h predict, with respect to MGN 08 and 
GH08, lower luminosities by a factor of 3-5. The theoreti- 
cal models are far from unique: other recent calculations of 
[Nell] line emission from a EUV-induced photoevaporative 
disk wind, that neglect X-ray irradiation, yield luminosities 



similar to thos e obtained by MG N 08 and GH 08 for their ref- 
erence models ( AlexandeJ2008 ). The observation of line shifts 
and broadenings, accessible through ground-based high resolu- 
tion spectrographs, may help to discriminate amon g the differ- 
ent proposed emission mechanisms (e.g. lHerczegetal.1120071: 
van Boekel etaDl2009l) . 

Detection of the [Ne II] 12.81/j,m line was first reported 
by Pascucci et al. ( 2007|) for four sta rs out of a sample of 6 
transition-disk systems. Lahuis et al. ('2007) detected the line 
in 15 more T-Tauri stars and EspaillaLet ak (20073 added three 
more detections. The relation of the Une strengths with X-ray 
luminosity and with other system parameters has remained un- 
clear: Pascucci et al- ( 2007.) report a correlation of the [Nell] 
line luminosity wi th Lx and an anticorr elation with mass ac- 
cretion rat^H Af ; 'EspaillatetaD (12007 '), complementing the 
Pascucci et al. (2007) data with their own, fail to confirm the 
correlation with L v and find a possibl e direct correlation with 
Af; the sample of ' Lahuis et al. f(l2007l) had only sparse X-ray 
data; in all cases the samples are small and inhomogeneous, 
comprising stars in different star-forming regions with differ- 
ent ages and distances. 

In this contribution we investigate the connections be- 
tween the neon fine structure line emission and the stel- 
lar/circumstellar properties with particular reference to the X- 
ray luminosity. We choose to focus on p Ophiuchi, one of 
the closest, youngest and most studied low-mass star form- 
ing regio ns in the solar neighborhood (for a recent review 
see WiUdng et al.l |2008). This is motivated by the fact that 
(i) the region has been extensively observed with Spitzer 
IRS, (ii) high quality X-ray data have recently been ob- 
tained by ourselves with t he Deep Rho Ophiuch i XMM-Newton 
Observation (DROXO, ISciortino et"al] l2006l) . allowing the 



cross-correlation of Spitzer sources with well-characterized X- 
ray emitters. Moreover, the young (<1 Myr) p Oph members 
have hard and luminous X-ray emission, characteristics that are 
expected to favor an observable disk response. YSOs in p Oph 
cover a range of evolutionary phases and include a significant 
number of Class I protostars. While on the one hand this fact 
makes our sample inhomogeneous, it also results in a better 
coverage of the star/disk/envelope parameter space. We here 
assume that p O ph is at a distance of 120pc, the most recent 
value derived by lLombardi et alJ (l2008h . 

This paper is organized as follows: we first introduce, in 
§|2] the main X-ray and NIR data, as well as additional data 
both original and from the literature; in §[3] we derive the quan- 
tities used for the subsequent analysis. We then compare the 
theoretical prediction for the X-ray ionization proxies with the 
observations (§21) and look for physically meaningful correla- 
tions between these and other stellar/circumstellar properties. 
Section|5]summarizes our results and presents our conclusions. 
An Appendix describes our method to characterize the YSOs in 
our sample by fitting theoretical models to their spectral energy 
distributions. 



' N ote, however, that the four [Nell] detections of IPascucci et alj 
span very small ranges of [Ne II] line and X-ray luminosities, 
both ~0.2 dex, comparable with uncertainties. The same applies to M, 
for which the range is 0.5 dex. 
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2. DROXO/Sp/fzer sample 

In order to correlate the [Ne II] and [Ne III] line strengths with 
the stellar X-ray emission and with the properties of the cir- 
cumstellar material, we decided to focus on a physically ho- 
mogeneous and well characterized sample of YSOs. Our sam- 
ple includes YSOs that: i) belong to the p Ophiuchi star form- 
ing region, and are therefore both young and relatively coeval, 
ii) have been observed with the Spitzer/IRS in high resolution 
mode, iii) fall in the field of view of DROXO, and have there- 
fore well characterized X-ray emission. 



2.1. Sp\tzer/IRS data 

We started by searching the Spitzer archive for IRS obser- 
vations of p Ophiuchi members in the XMM-Newton field 
(cf. i)2.2b performed with the short-high module (SH: A = 
9.9 - 16.9 Aim; R -600; sHt s ize = 4.7"xll.3'0 . Excluding 
GY65, which was identified bv iLuhman & Riekd fl999b as a 
background star, our sample consists of 28 YSOs. Table[T|lists 
these objects, the Spitzer program(s) under which they were 
observed, the total IRS (SH) exposure times and the details of 
the observing strategy: the number of exposures, the number of 
data collecting events (DCE) per exposure, and the integration 
time for each DCE. Note that four objects were targeted by two 
separate programs and have therefore been observed twice. 

We downloaded the short-wavelength, high-resolution 
Basic Calibrated Data (BCD) for the 28 stars in our sam- 
ple from the Spitzer archive. In order to produce final spec- 
tra we used the tools suggested on the Spitzer Science Center 
web page^. Specifically, after removing bad pixels with 
IRSCLEANw 1.9 we extracted the spectra of each DCE in the 
Spitzer IRS Custom Extraction (SPICE) v 2.0.4 environment. 
We then added up all the spectra from a given observation of 
a given target, from a minimum of two (the two nod positions) 
up to 72. This leaves us with 32 spectra (28 + 4 for the ob- 
jects that were observed twice). Finally, we used IRSFRINGE 
vl.l for the defringing. The background (sky) emission as a 
function of wavelength was estimated using spoi^ m steps 
of 0.5^m for A=10-19Aim and O.Ol^im for A=12.76-12.87/^m 
(the region of the [Ne II] line). These values, computed for the 
sky coordinates of the objects and the observation date, include 
the expected contributions from the Zodiacal light, the inter- 
stellar medium, and the cosmic infrared backgrouncfl They do 
not consider any eventual extended emission in the target prox- 
imity. Note, however, that this should not affect our main pur- 
pose, i.e. measuring the [Nell] and [Nelll] line fluxes, since 
emission from these lines is not expected from the cold molec- 
ular cloud in the absence of hot ionizing stars. Previous ob- 
servations of YSOs have moreover indicated that the emis- 
sion o f these lines is sp atially unresolved at the Spitzer reso- 
lution dLahuis et al.ll2007i) . Any continuum emission from the 
molecular cloud, if present, will thus be subtracted along with 
the stellar continuum when measuring line fluxes. However, 
multiple emission components eventually present within the 



Spitzer beam (4-5") will obviously be included in the [Nell] 
and [Nelll] fluxes, includin g e.g. those that may be associated 
with outflows as shown by Ivan Boekel et al.l (l2009b for the T 
Tauri system. 

2.2. X-ray data 

The Deep Rho Ophiuchi XMM-Newton Observation (DROXO) 
is the most sensitive X-ray exposure of the p Oph star form- 



ing region performed so far ( Sciortino et al. 2006h . It consists 



of an observation of Core F performe d with the XMM-Newton 



X-ray telescope dJansen et al.l 1200 Ih . The nominal pointing 
position was a2ooo = 16:27:16.5, (^2000 = —24:40:06.8. 
The observation, interrupted only by the constraints of the 
satellite orbit, was carried out in five subsequent revolutions 
(#0961... #0965), for a total exposure time of 515ksec. We use 
here data from the European Photon I maging Camera (EPIC; 



Striider et alJuOOlt [Turner et alj|2001h . consisting of three al- 
most co-pointed imaging detectors (MOSl, M0S2, and pn) 
sensitive to 0.3-lO.OKeV photons and with a combined field 
of view of '^0.2 square degrees. Source detection resulted in 
a list of 111 X-ray emitters, 60 of which are identified with 
a mid-infrared object detected bv lBontemps et alj (1200 lb with 
ISOCAM at 6.7 and/or 14.3/im. Details of the data reduction 
and general results from DROXO are found in Pillitteri et al. 
(2009, submitted). The present study is limited to the 28 YSOs 
with Spitzer/IRS coverage (see ! 12.1b . 

Twenty-two of our 28 YSOs are positionally matched with 
a DROXO source using a 5" identification radius. All the 
identifications are unambiguous. Cols. 5-7 of Table |4] list the 
DROXO source numbers from Pillitteri et al. (2009) and, for 
all 28 objects, off-axis angles and effective exposure times at 
the YSO position for the three EPIC detector^ Six YSOs are 
not detected in DROXO (see ji l3.1.2b . In three cases we used 
Chandra ACIS data from the literature. In the remaining three 
cases we have computed upper limits for the count rate as de- 
scribed in Pillitteri et al. (2009). Chandra ACIS data were also 
used for one of the DROXO-detected sources, IRS42/GY252. 
In the DROXO data the photon extraction region for this object 
is contaminated by a nearby bright source. The higher spatial 
resolution Chandra data is not affected by this problem. 

2.3. Ancillary data and SED fits 

We have collected additional data for our targets from which 
we derive relevant physical parameters. A summary of the re- 
sults is g iven in Table|2l Col. 2 h sts the ISOCAM source num- 
ber from lBontemps et aL ( 2001 ): col. 3 the YSO class derived 
both from the spectra l slope between 2 pm and 14 pm (reported 
from Bontemps et al.,.20 01) and from our own model fitting of 
the SED (see below and Appendix|A]l. As detailed below, the 
two classifications agree for 70% of the sources. 



http://ssc.spitzer.caltech.edu/postbcd/ irs-reduction.htmil 
http://ssc.spitzer.caltech.edu/propkit/spot 
http : //ssc . spitzer. cal tech. edu/ documents/background 



The off-axis angle gives an indication of the quality of the point 
spread function, which is sharpest on the optical axis; the effective ex- 
posure times are normalized values taking into account the vignetting 
of the optical system and the detector efficiency at the source position. 
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Table 1. Observation parameters for the SpitzerllKS and DROXO data. 



Object 


Prog.Id. 


Exp.Time 




DROXO* 


Off ax 


Exp.T 






[s] 


[s] 




['] 


[Ks] 


DoAr25/GY17 


179 


975.2 


4x2x121.9 


3 


12.4 


140/130/165 


IRS14/GY54 


179 


62.9 


2x1x31.5 




13.7 


127/ -/ - 


WL12/GYI1I 


172 


62.9 


2x1x31.5 


8 


9.1 


195/179/238 


WL22/GYI74 


93 


2265.0 


12x6x31.5 


27 


6.5 


251/234/285 


WL16/GYI82 


93 


251.7 


4x2x31.5 




4.2 


308/285/400 


WL17/GY205 


2 


12.6 


2x1x6.3 


34 


2.9 


326/317/440 


WL10/GY21I 


3303 


2265.0 


12x6x31.5 


35 


6.2 


263/ 17/337 


EL29/GY2I4 


93+2 


125.8 


6x3x6.3+2xlx 6.3 


38 


3.2 


309/309/434 


GY224 


172 


251.7 


4x2x31.5 


39 


1.4 


357/368/490 


WL19/GY227 


172 


251.7 


4x2x31.5 


40 


2.0 


359/341/469 


WL1I/GY229 


3303 


15414.1 


8x4x481.7 




5.4 


283/ 19/376 


WL20/GY240 


172 


251.7 


4x2x31.5 


46 


1.4 


376/357/486 


IRS37/GY244 


172 


62.9 


2x1x31.5 




11.2 


-/ 11/55 


WL5/GY246 


3303 


8776.6 


12x6x121.9 


49 


11.2 


-/1 1/206 


IRS42/GY252 


172+2 


75.5 


2xlx31.5+2xlx 6.3 


54 


1.9 


415/401/495 


GY253 


3303 


15414.1 


8x4x481.7 


55 


3.7 


367/390/478 


WL6/GY254 


172 


62.9 


2x1x31.5 


56 


10.3 


183/ 12/247 


CRBR85 


172+2 


306.7 


2x2x121.9+2x1x31.5 




2.0 


415/401/511 


IRS43/GY265 


2 


12.6 


2x1x6.3 


62 


2.5 


412/396/504 


IRS44/GY269 


2 


12.6 


2x1x6.3 


64 


2.6 


395/377/494 


IRS45/GY273 


179 


251.7 


4x2x31.5 


65 


13.0 


-/ -/1 84 


IRS46/GY274 


172 


251.7 


4x2x31.5 


67 


3.0 


388/374/487 


IRS47/GY279 


172 


62.9 


2x1x31.5 


68 


12.8 


-/ 10/196 


GY289 


3303 


15414.1 


8x4x481.7 


75 


7.6 


256/ 17/346 


GY29I 


3303 


3900.6 


8x4x121.9 


76 


8.4 


237/ 16/323 


IRS48/GY304 


2 


12.6 


2x1x6.3 




10.6 


196/ 13/263 


IRS5I/GY315 


172+2 


264.3 


4x2x31.5+2xlx 6.3 


87 


6.1 


345/356/432 


IRS54/GY378 


2 


12.6 


2x1x6.3 


97 


11.6 


187/ 13/257 



Column description. (1): Object names. (2) Spitzer programs from which the IRS data was taken. The identification numbers correspond to 
the following programs: 2=Spectroscopy of protostellar, protoplanetary and debris disks (P.I. : J.R. Houck); 93=Sun'ey of PAH Emission, 10- 
19.5pm (RI.: D. Cruikshank); il2l\19=From Molecular Cores to Planet-Forming Disks, (RI.: N. Evans); 3303=The Evolution ofAstrophysical 
Ices: The Carbon Dioxide Diagnostic (P.I.: D. Whittet). (3) Total IRS exposure time accumulated for each object in the short-high configuration. 
(4) Observing strategy, i.e. the number of exposures, ncxp, the number of data collecting events (DCE) per exposure, ndcc, and the exposure 
time of each DCE, Tdcc- For targets observed by multiple programs these figures are reported for each program. (5) DROXO source number 
from Pillitteri et al. (2009) for objects detected in the DROXO data. (6) Distance, in arcmin, from the XMM-Newton optical axis during the 
DROXO exposure. (7) Exposure times for the three EPIC detectors (PN/MOS1/MOS2); missing values indicate that the source was outside the 
detector FOV. 



Stellar parameters for Class II and III sources (according to 
our own classification based on SED fits) were estimated from 
the near-IR (2MASS) photome try. The procedure we h ave used 
follows closely that adopted by Bontemps et al. I (1200 1") and im- 
proved by Nf^Jg^Ud^ (2006). We assume that the J-band emis- 
sion from these sources is dominated by the stellar photosphere 
and that it is only marginally contaminated by the emission 
from circumstellar material and that the IR colors of Class II 
sources can be approximated b y emission from a p assive cir- 
cumstellar disk as described by Meyer et al.l (119971) . These as- 
sumptions obviously do not apply to Class I sources and for 
this reason we do not derive photospheric parameters for these 
sources. 

Our procedure starts with the dereddening of each object in 
the J — iJ vs. H — K diag ram to the locus of cTTS . As opposed 
to the procedu re used by 'N atta et al ( 2006 ^ we have used the 
Cardem et al.l (11989) extinction law. Two sources have colors 
slightly bluer than those of reddened main-sequence stars, pre- 



sumably due to photometric uncertainties. Dereddening these 
sources extrapolating the colors of Class II and III sources 
would produce an overestimation of the extinction. To mini- 
mize this effect, we have dereddened these sources to J — H = 
0.578. 

The values we derive for the J-band ext inction A j (col. 4 o f 
Table |2]l are very similar to the numbers in Natta et al. I (l2006h . 
The one exception is WL 16 for which our procedure produces 
a significantly higher extinction. 

Bolometric luminosities (col. 5 of Table |2]l were estimated 
from the dereddened J band magnitud es and the bolomet- 
ric correction used by iNatta et al.l (l2006h : logLboi = 1-24 + 
1.1 log Lj. Stellar masses and effective temperatures (col. 6 
and 7) were obtained fr om Lhoi assuming that s t ars lie on the 
0.5 Myr isochrone of the lO'Antona & Mazzitellil (Il997h evolu- 
tionary tracks. 

As part of the DROXO program, we have obtained comple- 
mentary IR spectroscopy at the VLT using the ISAAC instru- 
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ment a nd the same observing modes described in iNatta et al 
( l2006h . Low -resolution spectra (A/ A A 900) in the J or /iT 
band were obtained for 1 2 of the 13 YSOs in our sample that 
had not been observed by Natta et all (2006), the exception be- 
ing WL 19/GY277. These spectra comprise the Pa/3 and Br7 
lines that we use to measure accretion luminosities and mass 
accretion rates. For the reduction of the ISAAC data and the 
measurements of t he Pa/? or Br7 line we followed the proce- 
dures described by iNatta et al. I (l2006h . Accretion luminosities 
(or upper limits), b oth fro m the new near-lR spectra and from 
those o f iNatta etaLl(l2006h . were then estimated from empirical 
relations with the Pa/3 or Br7 luminosities (' Natta et alJl2004 ; 



Calvet et alJl2004 . They are hsted in col. 8 of Table IS Mass 



accretion rates (or upper limits), listed in col. 9 of the same 
table, were calculated from Lace and the photospheric param- 
eters derived from the near-IR photometry. They were there- 
fore computed only for YSOs for which these latter parameters 
are available, i.e. Class 11 and Class 111 objects with complete 
2MASS photometry. As a result the new near-lR spectra add 
only one accretion rate (for IRS 54/GY 378) and tw o upper lim- 



its (for GY 253 and IRS 5 1/GY 3 15) to the values in lNatta et al 
(l2006h . 



Given the fragmentary nature of the above described sys- 
tem parameters we have striven to obtain a more complete and 
homogeneous set of estimates by fitting the SEDs of our ob- 
jects with star/disk/envelope models. The details of the proce- 
dure, as well as the tests we have performed to ascertain its use- 
fulness, are described in the Appendix. Table [3] lists, for each 
source, the quality of the fit (the of the "best-fit" model) 
and the adopted values, with uncertainties for the following 
stellar and circumstellar parameters: extinction (the sum of in- 
terstellar and envelope extinction), stellar effective temperature 
and mass, disk mass, disk and envelope accretion rates. The 
last column indicates th e evolutionary stage as signed follow- 
ing the criteria given by iRobitaille et al.l ( l2007h . and reported 
in the Appendix ( ijA.2| i. based on the disk and envelope accre- 
tion rates and on the disk mass. These definitions are meant to 
reproduce in most cases the classical classification based on the 
SED slope (i.e. Class I, II, and III), which is often used to de- 
scribe the evolutionary status of YSOs. The stages, being based 
on physical quantities, have however the advantage of not de- 
pending on the inclination of the system with respect to the line 
of sight or on the effective temperature of the central object. 

Comparing the evolutionary stages from the SED fits with 
the IR classes derived from the ISOCAM photometry (Tab.|2] 
andO we find agreement in 19 out of 27 case^ (6, 11 and 2 
Class/Stage I, II, and III objects, respectively): 2 Class II ob- 
jects according to the ISOCAM classification are reclassified 
as Stage I and one as Stage III; 4 Class I and 1 Class III are re- 



As noted in the Appendix, gAT] the statistical significance of un- 
certainties is not easily assessed; the plausibility of the uncertainties 
on disk mass accretion rates is, however, indicated by a comparison 
with independent estimates from the literaure for a control sample of 
stars in the Taurus-Aurigae region (cf. Fig. lA.lb . 

' As described in the Appendix we rejected the result of the SED 
fits for one of our 28 YSOs, WL5/GY246. We consider it as a Class III 
object, i.e. the same Class given by the ISOCAM photometry, and 
derive its parameters from the spectral type and NIR photometry. 



classified as Stage II. In the following we will base our discus- 
sion on the evolutionary stages defined according to the results 
of the SED fits. However, in order to use a more familiar ter- 
minology, when referring to the 'stages', we will improperly 
adopt the term 'class'. 

3. Analysis 

We here describe the steps taken to derive X-ray and 
[Ne II]/[Ne III] luminosities from the XMM-Newton and the 
Spitzer IRS data, respectively. 

3. 1. X-ray luminosities 

We discuss separately the X-ray luminosities of the 21 YSOs 
for which an analysis of the X-ray spectra from the DROXO 
observation was possible (cf. ii l2.2l ). and those of the re- 
maining 7 objects for which we either make use of previous 
Chandra ACIS observations or we compute upper limits from 
the DROXO data. We will then discuss possible biases and un- 
certainties on the X-ray luminosities due to the high source ab- 
sorption and to their intrinsic variability. 

3.1 .1 . Spectral analysis of DROXO sources 

For the 21 YSOs with usable DROXO data the observed low- 
resolution X-ray spectra were fitted with simple emission mod- 
els convolved w ith the detector response using XSPEC v. 12.3.1 
(lArnaudlll9^ as described by Pilhtteri et al. (2009). We an- 
alyzed the time-averaged spectra accumulated during times of 
low background, i.e. excluding the intense background flares 
due to solar soft protons. This is the same time filter used by 
Pillitteri et al. (2009) for source detection, as it maximized the 
sensitivity to faint sources. It is not, however, the same filter 
used by Pillitteri et al. (2009) for spectral analysis. This latter 
differs from source to source and was devised to maximize the 
S/N by including times of high background when the source is 
bright enough to contribute positively to the S/N. Although the 
resulting spectra have higher S/N with respect to those based 
on the universal time filter we use here, the ensuing luminosi- 
ties are not suitable for our purpose as they do not scale linearly 
with the time-averaged luminosities. 

As done by Pillitteri et al. (2009) in many cases we fit- 
ted simultaneously data from all three EPIC instruments. In 
other cases the combined fits were statistically unsatisfactory 
because of cross-calibration issued and we excluded one or 
two of the detectors. The choice of detectors is the same as 
that of Pillitteri et al. (2009). In all cases, a model of isother- 
mal plasma emission (the APEC model in XSPEC) subject to 
photoelectric absorption (WARS) from material in the line of 
sight is found to be adeq uate. We adopted the plasma elemen- 
tal abundances derived bv lMaggio et al.l(l2007h for YSOs in the 



* One or both of the following: i) the source falls on a gap between 
the CCDs in one of the detectors, and we are unable to properly ac- 
count for the missing part of the PSF; ii) the source is intense and the 
statistical uncertainties per spectral bin are lower than the precision of 
the cross-calibration. 
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Table 2. Stellar/circumstellar data for the objects in our sample (see '>\23\ . 



Object 


ISO 


IR- 


SED^ 


A} 


1 L 

log — 


log Tcft 


M« 


log Lace 


log Mace 




Src. 


Class 


[mag] 


[K] 


[Mq] 


[Lq] 


[Mo/yr] 


DoAr25/GYI7 


38 


II 


-II 


0.7 


0.15 


3.63 


-0.27 


<-2.34 


<-9.23 


IRSI4/GY54 


47 


III? 


-III 


5.2 


-0.23 


3.55 


-0.57 


<-1.54 


<-8.15 


WLI2/GY1I1 


65 


I 


-I 


- 




- 


- 


- 


- 


WL22/GY174 


90 


II 


-II 


- 


~ 


- 


- 




- 


WLI6/GY182 


92 


II 


-II 


10.0 


2.26 


4.04 


0.56 


2.07 


-5.43 


WLI7/GY205 


103 


II 


-II 


11.3 


0.69 


3.68 


0.09 


- 


- 


WLI0/GY2I1 


105 


II 


-II 


4.5 


0.51 


3.67 


-0.04 


-0.89 


-7.92 


EL29/GY214 


108 


I 


-I 


- 


- 


- 


- 


- 


- 


GY224 


112 


II 


-II 


8.6 


0.54 


3.67 


-0.01 


- 


- 


WLI9/GY227 


114 


II 


-III 


16.3 


1.88 


3.90 


0.51 


- 


- 


WLI1/GY229 


115 


II 


-II 


4.3 


-0.94 


3.47 


-0.92 


-2.39 


-8.84 


WL20/GY240 


121 


II 


-I 


- 


- 


- 


- 


- 


- 


IRS37/GY244 


124 


II 


-I 


- 


- 


- 


- 


- 


- 


WL5/GY246 


125 


III 


-III 


16.8 


2.22 


4.02 


0.55 


- 


- 


IRS42/GY252 


132 


II 


-II 


7.5 


0.69 


3.68 


0.09 


<-I.14 


<-8.23 


GY253 


133 


III 


-III 


8.8 


0.36 


3.66 


-0.14 


<-I.53 


<-8.52 


WL6/GY254 


134 


I 


-II 


18.6 


2.43 


4.12 


0.59 


- 


- 


CRBR85 


137 


I 


-II 














IRS43/GY265 


141 


I 


-I 














IRS44/GY269 


143 


I 


-I 














IRS45/GY273 


144 


II 


-II 


6.6 


0.07 


3.62 


-0.32 


<-1.59 


<-8.45 


IRS46/GY274 


145 


I 


-I 














IRS47/GY279 


147 


II 


-II 


7.4 


0.63 


3.68 


0.05 


<-I.55 


<-8.61 


GY289 


152 


III 


-II 


7.3 


0.19 


3.64 


-0.24 


<-1.23 


<-8.14 


GY291 


154 


II 


-II 


7.4 


0.21 


3.64 


-0.23 


<-I.65 


<-8.58 


IRS48/GY304 


159 


I 


-I 














IRS51/GY3I5 


167 


I 


-II 


12.9 


2.29 


4.06 


0.57 


< 1.56 


<-5.95 


IRS54/GY378 


182 


I 


- II 


6.2 


0.35 


3.66 


-0.15 


0.02 


-6.95 



Notes: ^ IR class from the spectral slope between 2 ^m and 14 /xm jSontemps et alj200lh and SED Class from fitting of the SED (§|2j3] Table 

[3; 



Orion Nebula Cluster (ONC). In one case (EL29/GY214) the 
abundances of metals (all elements other than H and He) had to 
be incr eased by a factor of 3.5, with respect to the|Maggioetal 
(l2007h values, in order to obtain a reasonable fit. The spectra 
were fitted with a variety of initial parameters to avoid ending 
up in a local minimum of the space. The resulting fits are all 
statistically reasonable, with a mean Xred=l-1 ^'^'^ ^ maximum 
of 1.7. Results for the 21 usable DROXO sources (along with 
X-ray luminosities or upper limits for the other 7 YSOs; see be- 
low), are presented in TablelH cols. 6-9. For each source we in- 
dicate the detector(s) used for the spectral analysis, the Nh and 
kT values from the spectral fits, and the absorption coiTected 
X-ray luminosity in the 0.3 — 10 keV band. Statistical 90% con- 
fidence intervals for these quantities are also given. For Nh and 
kT these were obtained within XSPEC with the ERROR com- 
mand, while for Lx they were propagated from those on the 
plasma emission measures. 

3.1.2. Other objects 

We then estimated fluxes, or upper limits, for the 7 YSOs 
without a usable DROXO detection, marked in Table |4] by 
footnotes in the 'Instr' column indicating the source of the 
quoted iVn, kT, and ix values. IRS42/GY252 was detected 



in DROXO but, being contaminated by the wings of another 
bright source, we pr efer to use the luminosity obtained by 



Imanishi et alj (1200 Ih from the analysis of a Chandra ACIS 
spectrum, corrected for the different distance assumptions and 
choice of energy ban ds. IRS37/GY244 was also detected by 



Imanishi et alJ (1200 Ih but not in the DROXO data, possibly be- 



cause it lies close to the edge of the EPIC field and in the PSF- 
wi ngs of a brighter X-ra y source. The X-ray luminosity given 
bv llmanishi et alj ( l200lb for IRS37/GY244 is, however, based 
on a poorly constrained spectral fit and we, therefore, decided 
to estimate Lx from the ACIS count rate and a suitable conver- 
sion factor (see below). The same method was adopted to es- 
timate the Lx of IRS14/GY54 and WL16/GY182. T hese were 
detected as faint X-ray sources bv iFlaccomio et al. I (I2OO3.) in 



a re-analysis of the Chandra data but no spectral analysis is 
possible due to the low photon statistics. The remaining three 
objects, IRS48/GY304, WL11/GY229, and CRBR85, are not 
detected in any X-ray dataset and we estimated upper limits to 
their Lx from the upper limit to their XMM-Newton (DROXO) 
count-rate. 

Count-rate to luminosity conversion factors were thus em- 
ployied for six objects for which no reliable spectral anla- 
ysis was possible, i.e. three Chandra ACIS detections and 
three XMM-Newton upper limits. The conversion factors 
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Table 3. Parameters from SED fits(cf. §|A3. 



7 



Name 



2 

Xbest 



Av 
[mag] 



log Toff 
[K] 



log Mdisk 

[Mq] 



log Mdisk 
[Meyr-i] 



log Monv Stage/Class 
[Mej/r-i] 



DoAr25/GY17 


1 


01 


1.9^:76 


3.51^ 


0.33il:?g 


-2.0:^:? 




-7.6:;:^ 




—7.4-^'^ 


n 


IRS14/GY54 





11 


18.1?: 


3.5l:i 


0.35g:i 


-6.5ii:^ 


< 


-9.0 




— 


m 


WL12/GY111 


5 


46 


56.!i: 


3.4l:^ 


0.113:1? 


— 2.2_3 3 




""3-6.4 






I 


WL22/GY174 


31 


65 


/>Q 1.6e+03 
O'J-62. 


3.7l:^ 


l-5o:99 






-6.2i?:J 




— 8.3_8:3 


n 


WL16/GY182 


9 


53 


31.1: 


3.7|:? 


3.82:9 


_q 1 -1-2 




-8.0ll^ 




-7.8-''-^ 


n 


WL17/GY205 


2 


03 


42.1^: 


3.7l:^ 


1 7-^9 


-i.5il:^ 




-6.31^:^ 






n 


WL10/GY211 





51 


16.11: 


3.6l:^ 


l-Oo:42 


-3.2:^:^ 


< 


-8.3 




-8.2-''-^ 


n 


EL29/GY214 


6 


33 


42.i: 


3.83:7 


4.9i:^ 


-1-7:3:? 




-7.0:^:^ 




_c; 9-4.2 
•^•^-6.4 


I 


GY224 





33 


36.11: 


3.7l:? 


2.2?:i 


-3.7:^:? 


< 


-8.2 


< 


-7.9 


n 


WL19/GY227 


2 


70 


53.ii: 


3.83:8 






< 


-9.0 







m 


WL11/GY229 





04 


18.11: 


3.5l:i 


o.i5[]:ll 


-4.2_^:!^ 


< 


-8.9 




-8.7-'^-* 


n 


WL20/GY240 





81 


24.if: 


3.5l:i 


0.3Q[]:f? 


-2.7:y 




-7.9:|^ 




_c c-4.8 
"^■^-6.0 


I 


IRS37/GY244 


2 


02 


42.i: 


3.6l:^ 


0.63j:L 


— 2.2_3 7 




-8.i_g:g 




-5.91^:1 


I 


WL5/GY246t 




— 


52.44] 


Q q3.8 
•J-03.8 


3.o!:l 







— 







m 


IRS42/GY252 


1 


87 


30.i: 


4.it:J 


3.01:1 


-3-6:^:^ 


< 


-7.5 







n 


GY253 


5 


04 


29.28: 


3.6l:l 


o.7ig:l§ 


-8.0:^:^ 


< 


-9.0 






m 


WL6/GY254 


2 


10 


53.iE: 


4-^4.0 


Q q3. 5 
•J-J2.8 


0. / _4,7 


< 


-7.9 






n 


CRBR85 





25 


'oi-m. 


3.7l:^ 


9 n3 o 

^•Uo.79 


q 9-2.0 

■5-^-4.7 


< 


-7.9 




-6.51^:^ 


n 


IRS43/GY265 


1 


55 




3.6l:^ 


2.of:t 


-1 5-^'^ 




_7 1 -6.3 
' --^-8.5 




-5.ozt:t 


I 


IRS44/GY269 


2 


22 


57.1?: 


3.5l:i 


o.3i[]:i 


-2.21^:^ 




-6 1-5** 






I 


IRS45/GY273 





81 


29.il: 


q 74.0 

J- ' 3.6 


1-^0.97 


-4.91^:1 


< 


-9.0 


< 


-7.2 


n 


IRS46/GY274 


1 


23 


09 34. 


3.6l:^ 


0.6ll:le 


-2 l"^ "^ 




r, Q-6.9 

-'■0-9.2 




-5.8-?-° 


I 


IRS47/GY279 


4 


86 


3i.i: 


3.7l:? 


2.7'tt 


-1-9:^:^ 




-8.41^:^ 




-6.41?;? 


n 


GY289 


2 


32 


ryn 25. 

^'^•22. 


3.6l:l 


0.47j:lo 


-5.8:^:^ 


< 


-9.0 


< 


-7.0 


n 


GY291 





19 




3.6l:^ 


o.5oJ:l3 


-4.5:^:^ 


< 


-8.6 


< 


-7.5 


n 


IRS48/GY304 


4 


73 


26.11: 


3.6l:^ 


0.94§:i6 


__9 n-15 










I 


IRS51/GY315 


12 


91 


30.1?: 


3.7l:^ 


9-12.5 


-1 5-^'3 




-6.7:?:^ 




-8.1-'^-^ 


n 


IRS54/GY378 


1 


17 


26.11: 


3.7l:^ 


9 r,2.6 
^•^0.63 


15 

— 1.9_2 2 




-7.1:?:? 




7 9 — 6.5 
'■^-8.4 


n 



f : Parameters not derived from SED fits, but from the spectral type and NIR photometric (see Appendix). 



were cornputed with the Portable, Interactive Multi-Mission 
Simulatorij^PlMMS) assuming an isothermal plasma emission. 
This requires the assumption of a plasma temperature, kT, 
and, more crucially, of an absorption column density, A'h- 
For the temperature we took kT = 3.4 keV, the median of 
the kT values obtained from the X-ray spectral fits of the 
DROXO detections in our sample. Absorption estimates were 
derived from the Aj values in Table |2l when available, and 
from the Ay values in Table[3] in the remaining ca s es. A ^ 
and Ay were converted to A^h following I Vuong et al.l (120031) : 
Nh = 5.6 X 1021 Aj = 1.6 X lO^iAv cm'^. To the three Lx 
estimates from the Chandra ACIS count-rates we assign a 90% 
uncertainty of 0.5 dex. 

3.1.3. Biases and uncertainties 

Given the high absorption to which p Ophiuchi members are 
subject, we may wonder whether some or all of our X-ray lu- 
minosities are biased by the fact that low temperature emission 
components may be fully absorbed and therefore unaccounted 
for in the spectral fits. A very soft X-ray emission li ke that 
of the evolved CTTS TW Hya (e.g. IStelzer & Schmittl l2004). 
kT=Q.2-Q.3ke'V, would indeed have remained undetected in 



p Ophiuchi, as for a typical A^h = 4 x lO'^^cm^^ the ob- 
served flux in the XMM-Newton band is reduced by a factor 
~ 5 X 10'*, with respect to the unabsorbed case. The X-ray 
spectrum of TW Hya is, however, quite peculiar among YSOs. 
In the ~lMyr old ONC, for example, the Chandra Prior i 



http://heasarc.gsfc.nasa.govArools/w3pimms.html| 



Ultradeep Project (COUP) observation (iGetman et al.l l2005h 
indicates that, based on a sample of ~100 members subject 
to little absorption (A^h < 10^* cm^^) and whose X-ray spec- 
tra are well fit by 2-T models, the high-temperature component 
dominates the emission in most cases (80%) and, indeed, the 
mean kT (weighted by the emission measures of the two com- 
ponents) is >1.07keV in 95% of the cases. We therefore argue 
that, if low temperature components similar to those observed 
in the ONC should remain indeed unobserved with our data, 
the resulting underestimation of the X-ray luminosities would 
typically be less than a factor of two. 

Another source of uncertainty on the X-ray luminosities 
is their intrinsic time variability. While a full study of YSO 
variability in p Ophiuchi is beyond the scope of the present 
work, assessing its effect is important when we correlate Lx 
with other quantities observed non-simultaneously with the X- 
ray observation. Pillitteri et al. (2009) compare the average X- 
ray emission during the DROXO observation with that detected 
during previous Chandra and XMM-Newton observations of p 
Ophiuchi. The comparison indicates that the activity levels, av- 
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Table 4. Results from the analysis of the SpitzerllRS and DROXO datasets (cf. §|3TI]and ^ [T2] i. 



Name 


L[NcII] 


TTicont 
-T [NcII] 


L[NcIII] 


^ [NcIII] 


Instr. 


Nh 


kT 


Lx 




[10^** erg s-^ 


[Jy] 


[10^** erg s-^] 


[Jy] 




[lO^^cm^^] 


[keV] 


[10^** ergs^^] 


DoAr25/GY17 




0.26 


< 0.49 


0.28 


pn 


o.97S:;s^ 




180.64i?gi^ 


IRS14/GY54 


< 5.43 


0.28 


< 2.18 


0.12 


ACIS^ 






2.14?;g^ 


WL12/GY111 


43.29i?:i| 


39.20 


< 20.44 


28.55 


all 


14.62fi^* 


64.00 


12.14*«8f 


WL22/GYI74 


< 397.59 


73.63 


< 26.61 


11.97 


ml+pn 


i4.iiig:?2 


Qc:3.49 
/.OOi,72 


98.801^56" 


WL16/GY182 


< 80.07 


29.03 


< 36.24 


8.29 


ACIS^ 






^■iyo.69 


WL17/GY205 


< 23.40 


2.01 


< 17.84 


2.09 


all 


Q 7q5.91 
•J- ' -52.20 


3.43}°9?i 


7 qq17.53 
' ■003.97 


WL10/GY21I 


5.44|:i 


0.41 


< 2.79 


0.40 


pn 


2.80i:iSl 


4.85^:11 


39.66tg:i 


EL29/GY214 


< 78.98 


67.46 


< 40.82 


47.74 


ml 


6.37l:M 


4.23t:i? 


159.79?i?ig 


GY224 


< 5.51 


1.36 


< 4.33 


1.38 


ml+m2 


3.42tig 


8.62^306 


20.58f?if 


WL19/GY227 


< 15.17 


1.63 


< 9.32 


1.06 


all 


9.92ji# 


3.49ti 


7c; /I rl04.74 


WL11/GY229 


< 0.31 


0.03 


< 0.37 


0.04 


all^ 






< 2.69 


WL20/GY240 


19.80?g;gg 


1.74 


< 4.87 


3.09 


pn 


2.28i:?i 


2.48i:!f 


115.98lg|:l? 


IRS37/GY244 


20.49?ifi 


1.59 


< 4.72 


1.50 


ACIS^ 


5.8o!:gJ] 


> 1.70 


18.20i''7f 


WL5/GY246 


r ore. 45 
0-^04.05 


0.81 




0.38 


all 


6.62l:i!] 


4.05i:i? 


621.86|9e;g6 


IRS42/GY252 


< 7.73 


6.32 


< 6.89 


5.48 


ACIS"* 


3.90i:?J] 


1.30|| 


OQ 1 q72.44 
2».iSi2.02 


GY253 


< 0.63 


0.01 


< 0.43 


O.OI 


pn 


3.42i;i? 


2.622'44 


1 00 Q0I35.85 


WL6/GY254 


< 41.10 


28.77 


< 22.26 


19.70 


all 


r 7o7.36 
( 04.22 


S.OOfst^ 


40.2li?;gg 


CRBR85 


< 10.72 


4.98 


< 7.31 


4.21 


air 






< 7.41 


IRS43/GY265 


261.59ii:i 


12.30 


< 35.94 


12.44 


pn 


4.57t:IE 


3.ooiii 


077 nQ296.42 
Z i 1 .U0258.23 


IRS44/GY269 


90.62^^81' 


62.53 


< 87.12 


69.45 


pn 


7.i7^:if 


Q -71 4.04 
■J- ' I3.42 


242.14if^:i 


IRS45/GY273 


9.99jiif 


1.22 


< 2.82 


1.15 


all 


i.ooJ:ii 


A 77O.OO 
4- ' ' 1.59 




IRS46/GY274 


< 13.05 


6.63 


< 11.69 


5.44 


pn 


8.43i^63" 


4.78i:i 


14.72i^7f 


IRS47/GY279 


10.29^^61' 


4.74 


< 5.73 


3.81 


all 


i.95?:iJ] 


2.12?;i 


1 a OQ39.40 


GY289 


< 0.49 


0.05 


< 0.42 


0.06 


all 


i.93?:ig 


3.09i:ii 


27.79i:!i 


GY291 


< 1.01 


0.19 


< 0.70 


0.16 


all 


Qq2.60 
/.J^2.04 


2.50i:?I 


34.87i:i§ 


IRS48/GY304 


< 37.55 


15.80 


< 18.56 


23.48 


all^ 






< 25.70 


IRS51/GY315 


< 19.38 


7.23 


< 13.53 


6.31 


ml 


3.46ili 


2.67i:i 


111.89^f8?° 


IRS54/GY378 


< 19.34 


2.76 


< 8.87 


3.21 


all 


20.05^''4f 


2.32 


45.21J:'o'o'-" 



Notes: IR and X-ray fluxes and luminosities are corrected for extinction. Col. 6 indicates the XMM-Newton/EPIC or Chandra (ACIS) detec- 
tor(s) whose data was used for the fitting of the X-ray spectra: mI=MOS 1 , m2=MOS2, pn=pn, all=MOS l-l-M0S2-Fpn. Notes indicate the origin 
of A'h, kT, and Lx for objects with no usable DROXO detection (^13.1.21: ^ Lx fiom A CIS count-rate; ^ Lx upper limit from DROXO data; 
^ A^H and kT from th e spectral analysis of Chandra ACIS data bv llmanishi et al.l j200lh . Lx from ACIS count-rate; ^ A'h, kT, and Lx from 
Ilmanishietid] ( l200ll. Lx corrected for the differences in the assumed distance and energy band). 



eraged over ^^1 day, the typical length of the previous obser- 
vations, usually vary by less than a factor of 2 (Icr) over the 
timescale of years. The variability within each X-ray obser- 
vation, i.e. on the timescale of hours, can however be much 
larger due to flares that can reach up to ~100 times the qui- 
escent X-ray luminosity. These large flares are however not 
frequent. For the YSOs in the ONC, for example, an analy- 
sis of the hghtcur ves i n the COUP data s et alo ng the lines of 



Wolk et alj (l2005h and ICaramazza et al 



20071) indicates that 
^ by a factor of 2 or 



the X-ray flux is above the quiescent leve 
more for 10-15% of the time, and by a factor of 5 or more for 
2-3% of the time. Making the simplifying assumption that the 
Spitzer IRS observations are much shorter than the timescale 
of the X-ray variability, we can take these fractions as the frac- 
tions of objects for which the Spitzer observations coincided 
with an X-ray emission level above the characteristic level by 
more than the specified factor. For the 28 objects in our sample 
this implies 3-4 objects with a difference in Lx of a factor of 
>2 and < 1 with a difference of a factor of >5. 



'° M ore precisely the "characteris tic level" as defined bv lWolk et al.l 
( l2005h and lCaramazza et all bOOTh 



3.2. [Ne II] and [Ne III] line luminosities 

The estimation of neon line luminosities is performed in two 
steps: the direct measurement of fluxes from the reduced IRS 
spectra and the correction for extinction. 



3.2.1. Fluxes 

We measured the [Nell] and [Nelll] line fluxes, Fu and 
Fill, by integrating the spectra in the A= 12.79-1 2. 83/zm and 
A=15. 53-15. 57/im intervals, respectively. The underlying con- 
tinuum was subtracted by fitting a polynomial to two intervals 
on the left and right of the lines: 12.71-12. 78/xm and 12.84- 
12.9lAim for [Nell] and 15.45-15.53Aim and 15.58-15.65/Ltm 
for [Nelll]. The degree of the polynomial ranged between 1 
and 3 and the fit was repeated after excluding datapoints that 
deviated more than 2(t from a first fit. The la uncertainties on 
the fluxes, 6F11 and 6F111, were then estimated by propagat- 
ing the uncertainties on the individual spectral bins, taken as 
the maximum between the formal uncertainties given by the 
reduction process and the Icr dispersion of the continuum fit. 
The Unes were considered detected if the signal to noise ratio, 
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6F/F, was >3. In the opposite case upper limits were com- 
puted as max(i^, 0.0)h-3 SF. As indicated in §|TT] four YSOs 
have spectra from two separate observations: since the Ne Hnes 
are not detected in any of the spectra of these four objects we 
report the most stringent of the upper hmits and the continuum 
fluxes of the relative spectra. 

The [Nell] Une was detected in 10 YSOs, i.e. - 36% of 
our sample, while the [Ne III] line is detected only in one star, 
WL5/GY246, interestingly the brightest in X-rays in our sam- 
ple ( il4.2l l and likely a Class III object (see below). Figure [T] 
shows the 10 detected [Nell] lines and the single detected 
[Nelll] line. Gaussian profiles centered at the nominal line 
wavelengths and with normalizations from the measured fluxes 
are superimposed on the observed spectra. 

One of our YSOs IRS51 /GY 315, was also included in the 
study of Lahuis et al. ( 2007 *) using c2d data and we here use 



the same spectrum. Our 3a upper limits on [Ne II] and [Ne III] 
fluxes are ~20 % and ^50% high er than implied by the la up- 



per limits of iLahuis et al.l (l2007h . We attribute the discrepancy 
to the differences between the two measuring procedures. 

3.2.2. Extinction correction and luminosities 

In order to correct the line (and continuum) fluxes for extinc- 
tion we have chosen, for each of our YSO, a best-guess extinc- 
tion (Aj) from the up-to three estimates at our disposal. The 
Aj values from Tab.|2] estimated from the 2MASS photometry, 
were adopted when available. We otherwise estimated from 
the iVn values giv en by the X-ray sp ectral fits (Tab.Hji, con- 
verted according to lVuong et aP (l2003h (Aj = 1.8 • lO"^^ iVn). 
Finally, in the absence of the previous two estimates, we com- 
puted Aj from the Ay values gi ven by the SEP model fit s and 
listed in Tab.H Aj = 0.282 Ay dnieke & Lebofskvlll985l) . We 
make an exception to this rule for WL12/GY1 11, for which the 
iVn value is very uncertain and we prefer to use the extinctions 
from the SED fit. Table |5]lists these three estimates, along with 
uncertainties for the latter two, and the adopted Aj value. For 
extinctions taken from Tab.|2]we adopted a 1 cr uncertainty of 1 
magnitude. 

Critical for the derivation of unabsorbed fluxes is the choice 
of extinction-law, i.e., in the cases of the [Nell] and [Nelll] 
lines, the two ratios A12.S1 /Aj and Ai^^^^/Aj. The extinc- 
tion law at these wavelengths, in between two strong sili- 
cates absorption features at 9.7/im and 18^m, is not well 
established and seems to depend significan tl y on the grai n 



characteristics ( Weingartner & Drainel 1200 It iDraind 120031) . 



Chapman et al.l (120091) have recently established that, for stars 
in the p Ophiuchi region with low absor ption (Ak < 0.5), the 
Ry = 3.1 extinction law computed by IWeingartner & Draine 



(120011) for a "standard" grain size distribution fits the measure- 
ments between 1.25/im to 24/im. The extinction law of highly 
absorbe d stars (Ak > 2), however, is better reproduced by the 
i?i/=5.5 iWeingartner & Draine ( 200ll) extinction law, implying 
grain growth in the dense parts of the cloud. Since, with one ex- 
ception, the stars in our sample are highly extincted we adopt 
the Rv=5.5 extinction law and specifically A12.S1/ Aj — 0.16 
and A15.55/AJ = 0.13. 



0.080 
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0.25 
= 0.24 

Li_ 
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0.54 
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WL12/GY111 
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WL20/GY240 
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0.44 
0.42 

0.40 J 
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0.36 

0.34 
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15.8 



Fig. 1. Spitzer/IRS spectra in the [Nell] 12.81/im region for 
the ten YSOs for which the line was detected (upper panels), 
and in the [Nelll] 15.55/xm region for the single detection of 
this line (lower panel). The smooth (red) lines are the sum of 
the continuum fits used to measure the lines and Gaussians (of 
fixed width and cetenred at the nominal line wavelengths) with 
integral equal to the line fluxes measured by summing over the 
line spectral bins (see text). Object names are given in each 
panel. 
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Table 5. Different estimates of the J-band extinction for the 
objects in our sample, and adopted value (cf. ii l3.2.2| i. 



Name 




A 


7 [mag] 




Lit. 


Nh 


SED 


Adopted 


DoAr25/GY17 


0.7 


i- '1.6 


o.5!J:^ 


0.7 


1RS14/GY54 


5.2 






5.2 


WL12/GY111 




26.1 12.7 


15.812:8 


15.8 


WL22/GY174 






17.9l?6/ 


25.2 


WL16/GY182 


10.0 




8-67.8 


10.0 


WL17/GY205 


11.3 




11.91!:? 


11.3 


WL10/GY211 


4.5 


5.0!:? 


4.5t:J 


4.5 


EL29/GY214 




ii.4lg:! 




11.4 


GY224 


8.6 


f. 1 8.0 
6.I4.7 


10.2^°g« 


8.6 


WL19/GY227 


16.3 


17 720.6 


i4.9it:| 


16.3 


WL11/GY229 


4.3 




5.2i:g 


4.3 


WL20/GY240 




4-^3.9 


f. 79-9 
6./6.0 


4.1 


1RS37/GY244 




10.4^5/ 


1 r 812.7 
ii-OlO.8 


10.4 


WL5/GY246 


16.8 




i4.5li:i 


16.8 


1RS42/GY252 


7.5 


7.o!:i 


8.3^^2^ 


7.5 


GY253 


8.8 


6.ii:? 


8.i?:i 


8.8 


WL6/GY254 


18.6 


10.2^3^1 


i4.9l!:| 


18.6 


CRBR85 






i9.oli:i 


19.0 


1RS43/GY265 




o8.5 


i3.4lt:t 


8.2 


1RS44/GY269 




i2.8li:i 


i6.iil:i 


12.8 


1RS45/GY273 


6.6 


i.8g:i 


8.3?:? 


6.6 


1RS46/GY274 




iJ-ilO.l 


q i9-6 
^•^8.4 


15.1 


1RS47/GY279 


7.4 


3.5i:i 


79.4 
o-'8.3 


7.4 


GY289 


7.3 


, c3.9 
J-J3.0 


6.66.1 


7.3 


GY291 


7.4 


494.7 


0-^6.6 


7.4 


1RS48/GY304 






7 49.9 
'•^5.0 


7.4 


1RS51/GY315 


12.9 


6.25.5 


8.67.6 


12.9 


1RS54/GY378 


6.2 


35.8ii;6 


'•^7.0 


6.2 



With Aj ranging from 0.7 to 25 mag, the resulting cor- 
rection factors for the [Nell] 12.81/im line fluxes range from 
~1.1 to ^39 (mean: 5.9). Note that the difference with the 
Rv = 3.1 extinction law is significant: had we adopted it 
(Ai2.8i/Aj = 0.097) the 12.81/xm correction factor would 
have ranged from 1.1 to 9.5 (mean 2.70 

[Nell] and [Nelll] line luminosities were finally derived 
from the measured fluxes, assuming a distance of 120 pc. 
Resulting line luminosities and upper limits for the whole sam- 
ple are listed in Table HI along with absorption-corrected con- 
tinuum flux densities at the nominal line wavelengths. The re- 
ported uncertainties reflect measurement errors as well as un- 
certainties in Aj, but neglect possible systematic uncertainties 
related to the extinction law. 

We conclude this section with a cautionary note: the 
adopted extinction values refer to the central objects. The lumi- 
nosity corrections are thus valid only if the bulk of the [Ne II] 
and [Ne III] lines originates in the proximity of the YSOs. This 
assumptions may be false for emission from shocks associated 
with jets and outflows. 
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' ' The correction factors computed from the two extinction laws dif- 
fer significantly only for high extinctions: for DoAr 25, A] — 0.7, the 
difference is only ~4%. 



Fig. 2. [Nell] 12.81/im line luminosity vs. the continuum flux 
density at the same wavelength. Both quantities are corrected 
for interstellar extinction. Class I, Class II, and Class III objects, 
according to the classification derived from the SED fits, are 
indicated with different symbols (see legend). 



4. Results 

4.1. The [Nell] line 

We detect [Nell] 12.8 /im line emission in 10 out of the 28 
YSOs observed with Spitzer/IRS within the DROXO field 
of view (cf. Fig. [T]i- In one case, WL5/GY246, we also de- 
tect the [Nelll] 15.5/im line. All the [Nell] detections in 
pOph are X-ray sources: 9 are DROXO sources and the tenth, 
IRS37/GY244, was d etected in an earlier Chandra observation 
(llmanishi et al.ll200ll see i) l3.1.2b . Conversely the line is not 
detected in any of the three X-ray undetected objects. 

We investigated possible relations between the [Ne II] line 
emission and other physical parameters of the systems. First, 
however, we discuss an important observational bias, namely 
the dependence of our line detection sensitivity on the contin- 
uum intensity. Figure |2] shows the relation between the [Nell] 
line luminosity and the continuum flux density at the same 
wavelength. Both quantities are corrected for interstellar ex- 
tinction and stars of different evolutionary classes are plotted 
with different symbols. The lower boundary of detections and 
upper limits clearly shows a positive correlation, most likely 
due to the expected anti-correlation between the detection sen- 
sitivity and the counting-statistic uncertainties (that increase 
with the continuum). The upper envelope however also shows 
a correlation, which is independent of this detection bias and 
whose physical origin is to be understood. Figure |2] also indi- 
cates that Class I objects have higher continuum flux densities 
and [Ne II] line luminosities than Class II and Class III objects. 
This is not immediately interpretable in terms of the X-ray ex- 
citation mechanism discussed in the introduction. 
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Figure [3] shows the scatter plot between the [Ne II] line 
luminosity and Lx in 

the 0.3 - lOkeV banjH We also plot 
the six T-Tauri stars in the lPascucci et alj (l2007f) sample (four 
[Ne II] detections an d two upper lim i ts) and two stars (CS Cha 
and TWHya) from lEspaillat et aL ( 2007 ). Also shown are 
the theoretical predictions by 1VIGN08, ca l culate d, as a func- 
tion of Lx, assuming the iD'Alessio et alJ (Il999h disk model 

(M* = O.5A/0, i?* = 2.oi?0 , = mom, M = IQ-^Mq 

yr~^), and the predictions of GH08 for their fiducial model 
(model "A": M, = I.OMq, R., ^ 2.61Rq , = 4278K, 
M = 3 - 10"*Mq yr"\ Lpuv = 10^^ '^ erg/s) and two vari- 
ations: model "B" (with 100 times lower dust opacity, taken 
to represent an evolved disk) and model "D" (with 10 times 
higher FUV flux). 

Three conclusions are apparent: i) no overall trend of in- 
creasing [Nell] line luminosity with Lx is apparent; ii) two 
of the 10 measured [Nell] luminosities, those of DoAr25 and 
WL5/GY246, as well as the 18 upper limits, are consistent 
with predictions by current models for X-ray ionization of the 
disk; iii) the remaining 8 measured [Ne II] luminosities are 1 
to 3 orders of magnitude brighter than predicted. With respect 
to this latter point it is important to note, however, that the au- 
thors of both calculations stress that their models refer to ob- 
jects with specific star and disk parameters and are moreover 
affected by several uncertainties, e.g. in the atomic physics, in 
the simplified disk models that do not include holes, gaps, and 
rims, and in the current lack of EUV photoexcitation. Our find- 
ings indeed confirm these suspicions, and indicate that physical 
parameters other than Lx are likely to be important in deter- 
mining the [Ne II] line luminosity. 

Rather than a connection with Lx, Fig. [3] indeed suggests 
that the [Ne II] flux might be related to the evolutionary state of 
the YSOs, Class I being the strongest and Class III the faintest 



10' 



emitters. The position of TW Hya and of the six iPascucci et al 



(l2007h ' transition disk' -systems, seems consistent with this in- 
terpretation. CS Cha, also believed to host a transition disk has 
hower a strong [Ne II] emission in line with that of most of the 
[Ne II] detections in our sample. Within our sample, the higher 
[Nell] luminosity of Class I objects with respect to Class II 
ones is confirmed with significances ranging from 99.8% to 
99.99% by the five two-population tests for censored data im- 
plemented in the ASURV package ( Feigelson & Nelsonlll985 : 
[isobe & Feigelson 1990) . 

The overafl lack of correlation with X-ray luminosity leads 
us to investigate possible correlations of the [Ne II] emission 
with other stellar and circumstellar parameters. Figure|4]shows 
the relations with disk mass accretion rate and with stellar 
mass, both estimated from the SED fits. Also shown are the 
MGN08 and the GH08 model predictions. The most funda- 
mental of stellar parameters, mass, does not seem to influ- 
ence the [Nell] luminosity. At any given mass. Class I ob- 
jects have significantly higher line luminosities with respect 
to Class II and Class III ones; this points toward a role of 



• SED Class I (8) 
SED Class II (16) 

A SED Class III (4) 
Pascucci et al. (2007) 
CS Cha / TW Hya 




Lx [erg/s] 

Fig. 3. [Nell] luminosity vs. Lx- Symbols as in Fig. |2] The 
1 a error bars for the [Ne II] luminosity account for uncertain- 
ties in the measurement errors and in the extinction correction. 
For the X-ray luminosities estimated from DROXO data errors 
come from the statistical uncertainty of the normalization of 
the best-fit emission models, while for the Chandra ACIS de- 
tections we assumed a 50% uncertainty. Also plotted here are 
the six T-Tauri stars in the fPascucci et alJ (l2007l) sample and 
two stars from 'Es paillat et al.1 (l2007l) (CS Cha and TWHya). 
The model predictions of MGN 08 are indicated by the diago- 
nal fine; those of GH08, for their models 'A, 'B', and 'D' by 
the three diamond symbols. 



parameters related to the YSO evolution, such as matter in- 
flows and outflows. The [Nell] luminosity seems indeed to 
to correlate with M. Statistical tests for censored data, the 
Generalized Kendall's t and the Spearman's p as implemented 
in the ASURV package, confirm the existenc e of a correla- 
tion with ^99.5% confidence. The stars from IPascucci et al 



' The choice of energy band is not particularly important: had we 
restricted Lx to E> 0.87 keV, i.e. to photons relevant for the K-shell 
ionization of Ne, the points in Fig.[3]would have shifted downward by 
0.05-0.15 dex depending on the plasma temperature. 



(l2007h and lEspaillat et al.1 (l2007l) . shown in Fig.|4]but not used 
for the correlation tests, appear compatible with our sample. 
The correlation with Al may also explain the correlation of the 
[Nell] luminosity with the continuum flux at 12.81 fim (Fig.|2]i 
as this latter correlates strongly with M (not shown, confidence 
--99.99%). The [Nell] - M correlation might also explain, at 
least in part, the difference in the [Ne II] luminosity among the 
different evolutionary classes as Class I objects have statisti- 
cally higher accretion rates. 

As for the discrepancy with the X-ray excitation models, we 
note that three out of five Class I [Ne II] detections have nom- 
inal accretion rates that are higher than those used as inputs 
for both of the models considered here. This might be the rea- 
son for their higher than predicted line luminosities. The other 
Class I [Ne II] detections, WL 20 and IRS 37, however, have M 
estimates that, although with large uncertainties, are similar to 
those assumed by the X-ray ionization models and still their 
line luminosity is much larger than predicted. It is also pos- 
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sible that the [Ne II] -M correlation simply results from Class I 
objects having brighter line emission due to a mechanism unre- 
lated to disk mass accretion. One good candidate might be the 
defining characteristics of Class I objects, i.e. high envelope ac- 
cretion rates and/or their associated outflows. Indeed the [Ne II] 
luminosity also shows a significant correlation, at the ^99.9% 
level, with the Mem values derived from the SED fits. 

4.2. The [Ne III] line 

A spectral feature at ~ 15.55 /im, likely associated with a 
[Nelll] transition, is detected in only one star, WL5. An al- 
ternative identification for the observed feature might be a 
water rotational tran sition at 15.57^ 1, as detected e.g. by 
Can- & Najital (l2008h and ISalyk et alJ (l2008h in three CTTs. 
The line observed on WL5 is however well centered at 
15.55 iim (cf. Fig. [U and the wavelength difference with the 
water line is significant: ^2 spectral bins or about the FWHM 
spectral resolution. Moreover the m any other water lines that 
are seen in the spe ctra published by lCarr & Najital (2008) and 
Salvk et al. (2008') are not visible in the part of the WL 5 spec- 
trum shown in Fig.[Tl with the exception of a likely H2O line 
at 15.67 /im. We are therefore confident in the identification of 
the fine with [Nelll]. 

WL5 is a Class III objec t and, according to its F7 spectral 
type (IGreene & Meveiill995h . one of the most massive/hottest 
objects of the sample of [Nell] detections (see the Appendix 
for a detailed discussion of its properties). For both Neon lines, 
WL 5 has the lowest observed (i.e. absorbed) continuum flux 
among the stars in which [Ne II] was detected, thus facilitating 
the detections of the lines. The luminosities of both Neon lines 
compare reasonably well with the prediction of MGN 08 for 
the Lx of the object (6.2 x 10^° erg s~^): the [Ne II] fine is only 
7% fainter than predicted (well within la) and the [Ne III] line 
is 70% brighter than predicted (within 2a). 

All the other stars with [Ne II] detection in our sample have 
[Ne III] upper limits that are significantly larger and therefore 
compatible with the predictions of MGN 08. If we assume that 
the line ratios predicted by MGN 08, rather than the luminosi- 
ties, are correct and use the measured [Ne II] line luminosities 
to predict [Ne III] luminosities, we conclude that, for 8 out of 9 
stars, our detection sensitivity for [Ne III] is too low by a factor 
2.4-8. For the remaining case, IRS43/GY265, the star in our 
sample with the brightest [Ne II] line, the measured upper limit 
is only 10% higher than the predicted [Nelll] flux. 



5. Summary and discussion 

We investigated the origin of the [Ne II] and [Ne III] fine struc- 
ture lines by studying a sample of 28 p Ophiuchi members in 
the field of view of the DROXO deep X-ray observation and 
with available Spitzer IRS data. The [Nell] 12.81 /im and the 
[Nelll] 15.55 /im lines were detected in ten and one YSOs, 
respectively; absorption corrected line luminosities and upper 
limits for non-detections were computed and compared with 
predictions of X-ray disk ionization models. Finally, we ex- 
plored empirical relations between [Ne II] line luminosity and 



stellar and circumstellar parameters estimated by fitting the 
SEDs of the objects with star/disk/envelope models. 

The luminosities of the 10 detected [Nell] lines are, for 
the most part, 1-3 dex higher than predicted by models of X- 
ray irradiated (and ionized) circumstellar disks. Moreover, the 
[Ne II] luminosities do not correlate with the X-ray luminosi- 
ties. We conclude that, if these lines are indeed produced by 
X-ray ionization, factors other than ix are also important for 
the line production. Published models might still be valid: since 
they assume given star and disk characteristics (or few vari- 
ations) it is possible that some of these assumptions are crit- 
ical and that they do not correspond to the characteristics of 
most of our stars. Other excitation mechanisms might, how- 
ever, turn out to be more important than X-rays, such as strong 
shocks resulting from the interaction of the ste llar wind and 



jets with circumstellar mate rial (Hartmann & R aYmondlll989 
Hollenbach & McKeelll989"l:lvan Boekel et aLlbOoik 



Interestingly, the [Nell] luminosities of two of the objects 
in our sample, DoAr25 and WL5/GY246, match the theo- 
retical prediction for X-ray irradiated disks remarkably well. 
DoAr25 is a Class II object with mass accretion rate similar 
to that assumed by the models we have used for comparison. 
WL 5 is, based on its SED between 1 and 8 /im, a Class III ob- 
ject, but we cannot exclude the presence of a gas disk or a dust 
disk with a large inner hole. It mi ght be similar to the four stars 
with transitional disks studied by iBrown et al.l (12007 ), for one 



of which, T Cha, the [Ne II] line was detected by Lahuis e t al 
(20071). WL5 is, moreover, together with Sz 102 (Lahuis e t al 
I2OO7I) . the second YSO for which a detection of the [Ne III] 
line has been reported. As for the [Ne II] line, the luminosity of 
the [Ne III] line of WL 5 is roughly consistent with theoretical 
predictions for the X-ray irradiation mechanism. Also largely 
consistent with this mechanism are the upper limits to the lu- 
minosities of the undetected [Nell] and [Nelll] lines. 

The [Nell] 12.81 /im line luminosity correlates with the 
continuum flux at the same wavelength. While the lower enve- 
lope of the relation may be explained with a sensitivity bias, the 
upper envelope likely reflects a correlation with some physical 
characteristics to which the continuum flux is related. Disk ac- 
cretion rate, which we have found to correlate with the contin- 
uum flux at 12.81 /im, is one candidate, and indeed the [Nell] 
luminosity correlates with it. A tentative physical explanation 
of the correlation might involve the increased EUV flux pro- 
duced in the accretion shock, provided that this latter is able to 
reach and significantly ionize neon atoms above the accretion 
disk. Alternatively, given the correlation generally found be- 
tween accretion and outflow rates, the correlation might result 
from the [Nell] emission being produced in outflow -related 
shocks as mentioned above. Finally, the statistical correlation 
might be unphysical, and simply driven by the higher accre- 
tion rates of Class I objects combined with their higher [Ne II] 
luminosities. 

Indeed, Class I objects, i.e. those with significant envelope 
accretion, have significantly higher [Nell] luminosities than 
Class II objects. We propose that the presence of a circumstel- 
lar envelope and/or envelope accretion and/or the strong associ- 
ated outflows, i.e. the defining characteristics of Class I objects, 
plays a role in determining the line emission. Larger YSO sam- 
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Fig. 4. [Nell] luminosity vs. Af* (left panel) and Afdisk (right panel). Symbols as in Fig. |2] and [3] but with size proportional 
to logLx, as exemplified in the legend. The three X-ray undetected objects are indicated by empty symbols. JVIodel results are 
plotted as in Fig. [3] 



pies and more sophisticated theoretical models are needed to 
pinpoint the production mechanism of these gas-tracing lines 
and to derive a consistent picture of the environment around 
YSOs at different evolutionary stages. 
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Appendix A: (Circum)Stellar parameters from 
SED fits 

In this appendix we describe how we constrained some stel- 
lar and circumstellar parameters of the objects in our sam- 
ple by compari n g thei r SEDs with the theoretical models of 



Robitaille et al.1 (12006b . These consist of a grid of 200,000 
model SEDs that include contributions from the central star, the 
circumstellar disk, and the envelope, parametrized with 14 pa- 
rameters. The models that best approximate the observed SEDs 
were found wit h the a id of the Web based tool presented by 



Robitaille et al.l (120071) . As stated by Robitaille et al.1 (l2007h 



and in accord with basic principles, this method does not al- 
low the simultaneous determination of all the 14 physical pa- 
rameters, since the SEDs are often defined by less than 14 in- 
dependent fluxes. However, depending on the available fluxes, 
some of the parameters can be constrained more narrowly than 
others. We are here interested, in particular, in obtaining the 
range of values compatible with the observed SEDs for: i) the 
extinction toward our objects, ii) their disk accretion rates. 

A. 1. The method and its validation 



Our procedure follows closely that o f i RobitaiUe elaP tOQlh : 
from the Web interface we obtain, for each object, a list of 
the 1000 models that best approximate the observed SEDs, i.e. 
those with the smallest x^- Our "best guess" parameter val- 
ues and associated confidence intervals are then derived by se- 
lecting a set of statistically reasonable models and computing 
the median and the ilcr quantiles of the parameter values for 
these models. The statistically reasonable models were defined 
as those with reduced < (Xbcst + 3)' where Xbost refers to 
the best fit model, or in cases this condition results in less than 
10 models, the 10 models with smallest x^- Note that, because 
the uncertainties on the observed SEDs are not well defined 
(see below), and the parameter space is sampled only discretely 
by the adopted grid of models, the statistical significance of the 
thus derived confidence intervals cannot be easily assessed. 

A similar methoc(3 was tested by ^Rob itaille et al. ( 2007 ) 
by considering a sample of Taurus-Auriga objects for which 
stellar and circumstellar parameters had been derived indepen- 
dently in the literature and comparing these parameters with 
those obtained from fitting the SEDs, defined from the optical 
to millimeter wavelengths. In the case of our heavily absorbed 
p Ophiuchi YSOs the SEDs lack, with the exception of one 
star, data in the optical bands, i.e. those more directly affected 
by the accretion-shock emission. In order to test our ability to 
constrain the accretion rates in the absence of optical infor- 
mati on, we repeated the SED fits of the Taurus-Auriga stars 
of iRobitaille et al.l (120071) . using the same datapoints to define 
the SEDs, and both including and excluding the optical magni- 
tudes. Th e results are shown in F ig. I A. II Panel a), analogous to 
Fig 2b in IRobitaille et al.l (120071) . compares the accretion rates 
derived from the SED fits, including optical data, with indepen- 



IRobitaille et al.l ( |2007|) took the confidence interval for each pa- 
rameter as the full range of values in the selected models, and the 
"best guess" values as those of the model with minimum x^- 
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Fig. A.l. Compari son between mass accre tion rates from the literature and those derived from SED fits for the sample of T-Tauri 
stars considered in Robitaille et al. I (l2006l) . SED fits and determination of parameter ranges were performed as for the p Ophiuchi 
objects discussed in this paper. Panel a) compares the literature data with results of SED fits using all the available photometry, 
including optical bands. Panel h) is analogous, but only photometry longward of 1/zm was used for the SED fits. Panel c) finally 
compares the results of SED fits with and without optical photometry. Reduced values and mean absolute distances from the 
bisector, both computed considering uncertainties on the abscissa only, are reported within each panel. 



dent values from the hterature. Panel h) compares the results of 
the SED fits without the optical magnitudes with the literature 
data. The agreement between the two quantities is acceptable 
and may actually be considered better than in the former panel: 
the reduced x^, computed from the identity relation consider- 
ing only uncertainties on MsED, 'i^ indeed reduced from 
to 1 .7. This can in part be attributed to the increased error bars; 
note, however, that the average of the unsigned differences, 
abs{MsED-MLit), is almost unchanged, 0.49 dex for panel a) 
and 0.48 dex for panel h). Panel c) compares the M from the 
SED fits with and without optical magnitudes, showing that the 
two sets of values agree within uncertainties. We conclude that 
the SEDs defined from IR to milUmeter wavelengths are indeed 
sensitive to the accretion rate, at least in the M range covered 
by the Taurus-Auriga sample: log A/=[-8.5,-6]. 

This is due to the effect of viscous heating affecting the disk 
thermal structure. To exemplify this effect we plot in Fig. lA.2l 
as a function of accretion rate, t he ratio between the IR AC 3 
band and the J-band flux, for the lRobitaille et al. I (l2006h mod- 
els for stars with mass between 0.7 and 1.3 M©, age between 
1 and 2 Myr (implying little or no circumstellar envelope), and 
low disk inclination with respect to the line of sight {i <60°). 
We plot with different symbols models with disk inner radii 
in different ranges, since the inner hole affects the flux at the 
IRAC 3 wavelength (5.8 jim). A relation between the two quan- 
tities is seen for models with moderate inner disk holes, appar- 
ently characterized by different regimes in three different M 
ranges: \og{M/MQ) < -11,-11 < log{M/MQ) < -9, 
and \og{M / Mq) > —9. The factor of ~2 scatter around this 
relation may likely be attributed to model variations within 
the specified parameter ranges and to the several other uncon- 
strained model parameters. Similar and even more pronounced 
trends are apparent in analogous plots using fluxes in longer 
wavelength IRAC and MIPS bands, with the expected differ- 
ence that at the longer wavelengths, emitted farther out in the 



disk, the size of the inner hole has a much smaller effect. 
The three regimes in Fig. lA.2l can be understood as follows: 
i) for large accretion rates, log(Af/M0) > —9, the flux in 
the IRAC band, originated in the inner disk (R<1 AU), is sig 



nificantly affected by viscous accretion dD'Alessio et al.|[T998 
.1999.) ; ii) for -11 < log(A'//Mc 



< 



-9 disk heating is dom- 
inated by the stellar photospheric emission and, consequently, 
no relation between the IRAC flux and M is observed; in) for 
log(A//Afo) < —11 we again observe a direct relation be- 
tween the IRAC 3 flux and M, which we attribute to the fact 



that these low accretion rates correspond, in the lRobitaille et al 



(2006) model grid, to very low disk masses (Mdisk ^ 10~^M © 
for the '^l solar mass stars plotted in Fig. lA.2b . Since, in the 
model grid, disk mass and accretion are dire ctly correlated an d 
such low mass disks are optically thin (Robitaille et al.ll2006l) . 
lower accretion rates imply lower disk mass and lower emis- 
sion in the IRAC band. The IRAC 3 flux vs. M correlation in 
this regime does not therefore imply that that the mid-lR SED 
carries direct information on disk accretion. 

As a result of this discussion, in the derivation of accretion 
rates for our p Ophiuchi sample from the SED fits, we decided 
not to use values below lO^^Moyr^^. In such cases we in- 
stead conservatively assigned upper limits to M equal to the 
maximum between lO^^M © yr^^ and the upper end of the M 
confidence interval (see above). 

A.2. The p Ophiuchi sample 

We collected photometric measurements and uncertainties 
(when available) for our p Ophiuchi sample from several 
sources: J, H, and Kg magnitudes (or upper limits) were taken 
for almost aU objects from 2 MAS^j; Spitzer IRAC (bands 1- 
4) and MIPS (bands 1 & 2) photometry was collected from 



The J-band flux ofWL5 / GY24 6 and the ff-band flux of CRBR85 
were taken from lAUen etal] j2002h (converted from the HST bands); 
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Fig. A.2. Scatter plot of the ratio between the flux in the IRAC 1 
band ove r that in J, as a function of disk accretion rate, accord- 
ing to the Robitaille et al. ( 20061) mode ls for a solar mass stars . 
Each point corresponds to one of the iRobitaille et al. I (l2006h 
models satisfying the following conditions: mass of the central 
object between 0.7 and 1.3 M©, age between 1 and 2Myr, and 
disk inclination with respect to the line of sight <60°. Different 
symbols indicate models with inner disk radius in five ranges 
as indicated in the legend. 



the c2d databas tevans et al.' '2003*); 1.2 mm fluxes were 
collected from Stanke et al. (2006) and 1.3 mm fluxes from 
Andre & Montmerlel(ll994fR . Optical U BVR photometry for 



one object with small absorption (DoAr25) was taken from 



Yakubovl (Il992h . Table lATTI list all the photometric flux den- 



sities collected from the literature. 

Finally we complement the photometric data with flux den- 
sities from the IRS spectra (cf. i) l2.1l ). We computed flux den- 
sities between 10 and 18ptm, at regular wavelength intervals 
spaced by 0.5/im. Each flux density was taken as the average 
of the spectral bins in 0.2/im intervals centered at the nominal 
wavelength. For the four stars with two IRS observations we 
have taken the average of the two spectra. (In three cases the 
wavelength-averaged fluxes differ by less than 0. 1 dex, while 
in one case, EL29/GY214, the difference is 0.4 dex. In all cases 



the J- band upper limit for CRBR85 was taken from ISrandner et al] 

" Photometry extracted from the final (November 2007) c2d data 
delivery, selected according to the following conditions on quality 
flags: 'detection quality flag' equal to 'A', 'B', or 'U'; 'image type' 
equal to '0' for M1PS2 and different from '-2' and '0' for IRAC and 
MIPSl; 'flux quality' flag equal to 'A', 'B', or empty. 

Seven total fluxes fr om spatially resolved maps (Tab. 2 in 
lAndre & Montmerl3ll994) . and 15 peak fluxes or upper limits (from 
Tab. 1 in the same work), converted to total flux with the factors sug- 
gested by the authors. 



we verified that the results of the model fits did not change ap- 
preciably choosing either of the two spectra). Table |A.2l lists the 
flux densities from the IRS spectra. As stated in §|2T|our sky 
subtraction procedure does not take into account diffuse nebu- 
lar emission. In order to assess the significance of diffuse emis- 
sion on the object flux densities, we have considered the IRS 
spectra of the 13 YSOs in our sample observed in the context 
of the Spitzer legacy pro gram From Molec ular Cores to Planet- 



Forming Disks ('c2d'. lEvans et al. II2OO3'). As with the entire 



c2d sample, the reduced/sky-subtracted IRS spectra have been 
analyzed (and made publicly available) by the c2d team, us- 
ing a sophisticated extraction and sky subtraction method based 



on th e modelling of the cross dispersion profiles (iLahuis et al 



20071). We have compared the flux densities derived from the 



c2d-reduced spectra with those derived from the same spectra 
reduced by us. We find the spectra to be similar, with both the 
maximum and the wavelength-averaged discrepancy decreas- 
ing with object intensity. The maximum discrepancy falls be- 
low 10% for the 9 YSOs with c2d-reduced spectra that have av- 
erage flux >0.5 Jy. Based on this comparison, and noting that 
the c2d objects are representative of our sample as for their po- 
sition with respect to nebulosity seen in IRAC and MIPS maps, 
we decided to use the IRS-derived fluxes for defining the SEDs 
of the 17 stars with average I RS flux >0.5 Jy. 

As suggested by Robitaille et al.l(l2007h . in order to account 
for systematic uncertainties, underestimation of the measure- 
ment errors, and intrinsic object variability in time, a lower 
limit of 25%, 10%, and 40% was imposed on the uncertainties 
of optical, NIR/MIR, and millimeter fluxes, respectively. 

Figure lAJj exemplifies the "fitting" procedure described in 
^ lA.ll for three of our YSOs. It shows the SEDs with the best 
fit models and the distributions of two fit parameters. Ay and 
Mdisk, both for the 1000 models with lowest and for the 
statistically reasonable ones (cf. lA.ll ). SEDs and best fit models 
for all the 28 YSOs in our sample are shown in Fig. lA.4l 

Following visual examination of the SED fits and of the 
distributions of model parameters used to define the confi- 
dence intervals, we have decided to modify the input data- 
points for two objects: for IRS45/G Y273 we have excl uded 
the 1.2 and 1.3 mm d atapoi nts from IStanke et al.l (l2006l) and 
lAndre & Montmerle ( 1994 ): in both cases including these 
points significantly worsened the quality of the fit and had a 
significant effect on the values of the parameters. The 1.2 mm 
flux is >20 times higher than the 1.3 mm flux (an upper limit) 
and probably refers to an extended source that includes our 
YSO. For GY289, a source with average IRS flux <0.5 Jy, we 
have decided to include the IRS datapoints because: i) they 
agree quite well with the MIPS fluxes at similar wavelengths, 
ii) the quality of the model fit is reasonable (x^cst 
Hi) the confidence intervals of the model parameters are nar- 
rower but compatible with those from the fit performed without 
these points. 

Finally for one object, WL5/GY246, we could not obtain 
a unique fit with the above procedure. The object was previ- 
ously classified a s a deeply absorbed Cl ass III star with an 
F7 spectral type (iGreene & Meyeiill995l) . and our SED was 
defined by J,H,K, Spitzer IRAC 1-4 and 1.2/1. 3mm fluxes. 
Fits both with and without the mm fluxes, Ukely contaminated 
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Fig.A.3. Examples of SED fits for three objects in our sample with [Nell] detections. From left to right: DoAr25/GY17, 
WL20/GY240, and IRS44/GY269. The first is classified as Stage/Class II, the o ther two as Stage/Class I. The upper row shows 
the SEDs and the best fit models as produced by the Web interface provided by Robitaille et alJ (2006^. For the datapoints, de- 
tections and upper limits are indicated by circles and triangles, respectively. The lower two rows represent distributions of two 
fit parameters. Ay and Mdisk- The empty histograms refer to the 1000 model fits with lowest ™d the green histograms to the 
statistically reasonable samples of models defined in ij lA.ll The solid and dashed vertical lines indicate the median and the la- 
dispersion for these latter samples. For the panels in the second row the symbols close to the upper axis indicate the Ay values 
inferred from the Aj in Tab le|2] (circles) and from the X-ray derived A^h in Table|4](squares). 



by ne arby sources (c.f. Andre & ]VIontmerlell994 ; Stanke et al 
2006h . consistently yield high envelope and/or disk accretion 
rates, typical of a Class I object, but having little effect on 
the NIR/MIR part of the SED due to the associated large in- 
ner disk radii. The NIR/MIR SED can however be fit equally 
well by pur ely photospheric 'Phoe nix' models as suggested 
by the same iRobitaille et alJ (l2007h web interface used to fit 
the star/disk/envelope models. We thus decided to assume that 
WL5/GY246 is a Class III object and to derive its extinction, 
effective temperature, and stellar mass using the J, H, and 
K photometry, the spectral ty pe, and the calibrations tabulated 
by iKenyon & HartmannI (119951) . Uncertainties were estimated 
from the assumed uncertainty on the spectral type, one sub- 
class, and the range of values obtained by estimating the ab- 
sorption from the J-H, H-K, and J-K colors. 

Table [3] introduced in the main text ( i]2.3b lists the out- 
come of the SED-fit process: the quality of the fit (the of 
the "best-fit" model), the object extinction (the sum of inter- 
stellar and envelope extinction), the stellar effective tempera- 



ture and mass, the disk mass, the disk and envelope accretion 
rates, the evolutionary Stage. Th is latter quantity was assigned 
following liobitaille et al.1 (l2007h . Stage I: Menv/M* > 10"^; 
Stage II: Monv/M* < 10"^ and Mdisk/M* > 10"*^; Stage III: 
Mcnv/M^ < 10"*^ and Mdisk/M* < 10"*^. As indicated in 
the main text, in order to use a more familiar designation to re- 
searchers in the field, we also refer to the 'Stages' as 'Classes'. 



Figures IA.5I and IA.6I compare the extinction values (Ay) 
and stellar Tcff obtained from the SED fits with the same pa- 
rameters listed in Table|2]for Class II and Class III stars. Given 
the considerable uncertainties of both determinations, the SED 
fits yield results s imilar to those obtained with the method of 
Natta et al. ( 2006h . A similar comparison with the accretion 



rates derived from the Pa/? and Br7 NIR line fluxes (in Table 
|2]), is less conclusive due to the large number of upper lim- 
its and to the large uncertainties that affect the spectroscopic 
measurements as well as the SED fits. Seven objects can be 
used for the comparison having accretion rate estimates or up- 
per hmits from both methods. For only two stars both meth- 
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Fig. A.5. Comparison o f the Ay values derived from fitting the 
SEDs with the Robitail le et al.l 12006) models, with values de- 
rived from 2MASS photometry (cf. Table|2]i. Objects of differ- 
ent SED Class are indicated by different symbols as shown in 
the legend. 



(12006 . see also j]2.3b is better suited for cool stars and likely 
to yield inaccurate r esults for WL 16 {T pf; ^ K). An inde- 
pendent estimate by^Naj ita et al. (1996) yielded an upper limit 
compatible with the SED value: M < 2 x 10" MqJT^'^. Other 
three stars have Al estimates from the SED fits and upper limits 
from Table|2] in two cases, IRS 51 and IRS 47, the confidence 
intervals from the SED fits are consistent with the upper lim- 
its; for DoAr25/GY17, the only star with optical magnitudes, 
the SED fit yields an accretion rate that is 1 .6 dex higher than 
the upper limit from the Pap line. Finally, for two stars, WL 10 
and WL 1 1 , the spectroscopic estimates are 0.4 dex and 0. 1 dex 
larger than the upper limits from the SED fits. The discrepancy 
is however reduced to 0.24 dex for WL 10 and disappears for 
WL 1 1 if the slightly larger M values from lNatta et al.1 (|2006|) 
are considered instead of those in Table |2] 

A. 3. Summary 

In this Appendix we have shown that the SED models of 
Robitaille et al. (l2006h . although undeniably approximate, can 



be useful to constrain parameters such as the line-of-sight ab- 
sorption and the disk accretion rate, even in the absence of op- 
tical photometry. Although resulting uncertainties in these pa- 
rameters are often large, the constraints are by and large com- 
patible with independent determinations obtained with more 
direct methods. 
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Fig. A.6. Same as Fig. lA.5l for the effective temperatures. 



ods yield estimates: those for IRS 54 are in good agreement; 
for WL 16 the spectroscopic estimate is 2.6 dex higher than the 
value from the SED fits, M 10~^ MQyr~^. The discrepancy 
is however reduce d to 1 .2 dex wh en comparing the result of the 
SED fit with the .Natta et alj ( 12006) value. Moreove r the deriva- 
tion of M from the Pap line with the method of iNatta et al 



Table A.l. Flux densities, in mJy, collected from the literature (cf. ^A.2i and used for the SED fits. 



Name 






J 




H 




Ks 


IRAC 


[1] 


IRAC [2] 


IRAC 


[3] 


IRAC 


[4] 


MIPS 


[1] 


MIPS [2] 


1.2mm 


1.3mm 


DoAr25/GY17 


279. 


00 


448. 


,00 


484, 


,00 


367, 


,00 


292.00 


299, 


,00 


258, 


,00 


399 


.00 




153 


.00 


280, 


.00 


IRS14/GY54 




1. 


,66 


12. 


10 


30, 


.20 


35, 


,10 


29.20 


23, 


,40 


14, 


,00 












< 30, 


.00 


WL12/GY111 




0, 


,32 


3, 


,21 


15, 


.20 


239, 


.00 


744.00 


1610 


.00 


2240, 


,00 






8120.00 


415, 


.00 


130 


.00 


WL22/GY174 


< 


0, 


,20 


< 0, 


,08 


0, 


.59 












3230, 


,00 












400, 


.00 


WL16/GY182 




3, 


,44 


65, 


,90 


397, 


.00 


1400, 


.00 


1970.00 


5030, 


.00 










29100.00 






< 10 


.00 


WL17/GY205 


< 


0, 


,19 


1, 


,96 


27, 


.30 


240, 


.00 


416.00 


553, 


.00 


695, 


,00 


2790, 


.00 


6070.00 


144, 


.00 


70, 


.00 


WL10/GY211 


15. 


,30 


86, 


,10 


181, 


.00 


259, 


.00 


310.00 


272, 


.00 


222, 


,00 


339, 


.00 


784.00 






< 30, 


.00 


EL29/GY214 




0, 


,31 


39, 


,00 


929, 


.00 








12800, 


.00 












316, 


.00 


300, 


.00 


GY224 


< 


0, 


,20 


7, 


,21 


55, 


.70 


203 


.00 




358, 


.00 


367, 


,00 


907, 


.00 


908.00 










WL19/GY227 


< 


0, 


,06 


0, 


,97 


25, 


.20 


215, 


.00 




406 


.00 


328, 


,00 


223, 


.00 








< 20, 


.00 


WL11/GY229 




0, 


,90 


5, 


,83 


17, 


.00 


33, 


.70 


35.80 


34, 


.20 


31, 


,00 


43, 


.00 








< 20, 


.00 


WL20/GY240 




4. 


,41 


32, 


,00 


97, 


.30 


127, 


.00 


143.00 


140 


.00 
















90 


.00 


IRS37/GY244 


< 


0, 


,21 


1, 


,75 


15, 


.50 


127, 


.00 


206.00 


286, 


.00 


268, 


,00 


780, 


.00 








300, 


.00 


WL5/GY246 




0. 


02 


1. 


37 


39. 


90 


209. 


.00 


297.00 


298. 


00 


163. 


00 








<375. 


00 


35. 


00 


IRS42/GY252 




1, 


,31 


32, 


,50 


270, 


.00 


1060 


.00 


1630.00 


2100 


.00 


2980, 


,00 


3450, 


.00 


2940.00 






< 35, 


.00 


GY253 




0, 


,20 


5, 


,39 


32, 


.70 


64, 


.70 


60.40 


53, 


.30 


30, 


,90 


3 


.18 












WL6/GY254 


< 


0, 


,06 


0, 


,72 


31 


.10 






925.00 


1440, 


.00 


1730, 


,00 


4360, 


.00 


5110.00 






< 75, 


.00 


CRBR85 


< 


0, 


,01 


0, 


,06 


2, 


.26 


54, 


.00 


131.00 


191 


.00 


198, 


,00 


1340, 


.00 




< 921 


.00 






IRS43/GY265 


< 


0, 


,06 


4, 


,01 


84, 


.30 


629, 


.00 


1240.00 


1790, 


.00 


2190, 


,00 








967, 


.00 


190 


.00 


IRS44/GY269 


< 


0. 


,38 


3, 


,46 


47, 


.00 


731, 


.00 


1830.00 














34700.00 


< 576, 


.00 


180 


.00 


IRS45/GY273 




0, 


,81 


12. 


,20 


60, 


.70 


187, 


.00 


272.00 


382, 


.00 


481, 


,00 


712, 


.00 




1058. 


.00 


< 40 


00 


IRS46/GY274 




0, 


,33 


9, 


,03 


77, 


.50 






271.00 


402, 


.00 


411, 


,00 












90 


.00 


IRS47/GY279 




1, 


,19 


25. 


,40 


164, 


.00 


740, 


.00 


1190.00 


1580, 


.00 


2040, 


,00 


1720, 


.00 




1336, 


.00 


< 20, 


.00 


GY289 




0. 


,56 


8, 


,24 


29, 


.10 


45, 


.60 


41.90 


36, 


.80 






56, 


.50 












GY291 




0. 


,53 


8, 


,18 


27, 


.50 


54 


.20 


69.50 


84, 


.80 


108, 


,00 


94, 


.50 












IRS48/GY304 


94, 


,60 


305, 


,00 


618, 


.00 


1410, 


.00 


1600.00 


4060 


.00 


6000, 


,00 












240, 


.00 


IRS51/GY315 




0, 


,24 


14, 


,40 


169, 


.00 


752, 


.00 


916.00 


1000 


.00 


1070, 


,00 


2730, 


.00 


3260.00 


124, 


.00 


165 


.00 


IRS54/GY378 




2, 


,14 


34, 


,30 


218, 


.00 






712.00 


931, 


.00 


1010, 


,00 


3560, 


.00 


6500.00 


340, 


.00 


< 120 


.00 



Note: Values in italic were not used for the SED fits. 



Table A.2. Flux densities, in Jy, obtained from the IRS spectra for the SED fits. 



ISl amp / A 1 /f m 1 


10.0 


10.5 


1 1.0 


1 1 5 


12 


12.5 


13.0 


13.5 


14.0 


14 5 


1 5 


15.5 


16.0 


16 5 


17 


17.5 


18.0 




23 


23 


24 


27 


25 


25 


25 


23 


24 


22 


24 


2^ 


28 


50 


33 


36 


37 


1RS14 


022 


022 


040 


094 


096 


11 


098 


081 


066 


061 


062 




07^ 


11 


11 


09'y 


067 


WL12 


93 


1 3 


1 7 


2 3 


3 


3 7 


4 1 


4 3 


4 5 


4 7 


4 1 


A 1 


A A 


4 2 


4 3 


A A 


A 5 


WL22 


16 


25 


62 


87 


1 


1 4 




I I 


96 


96 


44 


48 


96 


1 5 


1 3 


1 4 




WL16 


1.3 


1.5 


4.4 


5.1 


4.8 


6.3 


4.4 


3.8 


3.2 


2.9 


2.6 


2.4 


2.9 


4.4 


4.0 


3.9 


2.7 


WL17 


0.18 


0.21 


0.27 


0.31 


0.34 


0.37 


0.42 


0.48 


0.48 


0.52 


0.51 


0.54 


0.64 


0.64 


0.70 


0.74 


0.76 


WLIO 


0.17 


0.18 


0.19 


0.20 


0.21 


0.21 


027 


022 


022 


022 


025 


0.24 


025 


026 


027 


0.27 


02S 


EL29 


7.2 


9.5 


12. 


16. 


19. 


22. 


23. 


24. 


24. 


26. 


22. 


22. 


24. 


24. 


24. 


25. 


25. 


GY224 


0.27 


0.30 


0.34 


0.37 


0.38 


0.39 


0.40 


0.43 


0.45 


0.47 


0.48 


0.48 


0.53 


0.56 


0.60 


0.62 


0.63 


WL19 


0.053 


0.071 


0.077 


0.096 


0.13 


0./5 


0.15 


0.16 


0.17 


0.18 


0.16 


0.15 


0.16 


0.14 


0.14 


0.15 


014 


WLll 


0.012 


0.014 


0.011 


0.010 


0.014 


0.014 


0.015 


0.019 


0.017 


0.017 


0.022 


0.021 


0.027 


0.029 


0.034 


0.037 


0.033 


WL20 


0.34 


0.45 


0.59 


0.73 


0.83 


0.89 


1.0 


1.2 


1.4 


1.6 


1.6 


1.8 


2.3 


2.5 


2.8 


3.1 


3.4 


mS37 


0.13 


0.17 


0.22 


0.26 


029 


05J 


035 


057 


059 


0.43 


0.40 


0.42 


049 


055 


057 


065 


066 


WL5 


0.015 


0.021 


0.029 


0.037 


0052 


0.064 


0.062 


0066 


0065 


0.066 


0.047 


0.047 


0.068 


O079 


0.088 


0.11 


0.098 


IRS42 


1.8 


2.0 


2.3 


2.4 


2.1 


2.1 


2.1 


2.1 


2.1 


2.2 


2.1 


2.1 


2.4 


2.4 


2.4 


2.6 


2.7 


GY253 


9.4e-05 


0.00031 


-0.0015 


-0.0021 


0.0024 


0.00084 


0.0028 


0.0059 


0.0035 


0.0023 


0.0051 


0.0024 


0.0047 


0.0048 


0.0080 


0.0092 


0.0058 


WL6 


0.73 


0.92 


1.1 


1.4 


1.6 


1.8 


2.0 


2.1 


2.2 


l.A 


1.9 


2.0 


l.A 


2.5 


2.6 


1.1 


2.8 


CRBRS5 


0.066 


0.10 


0.14 


0.18 


0.23 


0.27 


0J7 


055 


0.39 


0.43 


055 


0.37 


0.49 


057 


OJJ 


0.62 


0.65 


IRS43 


0.87 


1.2 


1.7 


2.3 


2.9 


3.5 


4.0 


4.4 


4.9 


5.4 


4.0 


A.A 


5.6 


5.8 


6.3 


6.8 


7.3 


IRS44 


1.5 


2,4 


3.5 


5.3 


7.2 


8.7 


10. 


12. 


13. 


15. 


13. 


14. 


17. 


18. 


20. 


22. 


24. 


IRS45 


0.34 


0.40 


0.44 


0.48 


0.46 


0.48 


0.48 


0.50 


0.52 


0.53 


0.49 


0.52 


0.59 


0.57 


0.60 


0.64 


0.63 


IRS46 


0.31 


0.39 


0.48 


0.57 


0.66 


0.72 


0.76 


0.81 


0.83 


0.87 


0.85 


0.88 


0.93 


0.96 


1.0 


1.0 


1.1 


IRS47 


1.2 


1.3 


1.5 


1.6 


1.6 


1.6 


1.6 


1.6 


1.7 


1.7 


1.4 


1.5 


1.7 


1.6 


1.7 


1.7 


1.8 


GY289 


0.0039 


0.0065 


0.0080 


0.011 


0.013 


0.014 


0.015 


0.021 


0.018 


0.019 


0.024 


0025 


0.030 


O055 


O057 


0.041 


O056 


GY291 


0.055 


0.060 


0.063 


0.064 


0.062 


0.062 


0.063 


0.065 


0.063 


0.064 


0.064 


0.063 


0.075 


O0S5 


0.091 


0.095 


0.091 


IRS48 


l.A 


2.6 


4.3 


A.l 


4.5 


4.9 


4.6 


5.1 


5.8 


6.7 


8.0 


9.2 


12. 


14. 


16. 


18. 


20. 


IRS51 


1.2 


1.3 


1.3 


1.3 


1.3 


1.2 


1.2 


1.3 


1.3 


1.4 


1.1 


1.3 


1.7 


1.8 


2.0 


2.1 


2.2 


IRS54 


0.48 


0.59 


0.73 


0.89 


1.0 


1.1 


1.2 


1.2 


1.3 


1.4 


1.4 


1.5 


1.7 


1.7 


1.8 


1.9 


2.1 



Note: Values in italic were considered uncertain and were not used for the SED fits. 



